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Abstract 

We present an expression for the leading-color (planar) four-loop four-point amplitude of = 4 
super symmetric Yang-Mills theory in 4 — 2e dimensions, in terms of eight separate integrals. The 
expression is based on consistency of unitarity cuts and infrared divergences. We expand the 
integrals around e = 0, and obtain analytic expressions for the poles from through 
Wc give numerical results for the coefficients of the and poles. These results all match 
the known exponentiated structure of the infrared divergences, at four separate kinematic points. 
The value of the coefficient allows us to test a conjecture of Eden and Staudacher for the 
four-loop cusp (soft) anomalous dimension. We find that the conjecture is incorrect, although our 
numerical results suggest that a simple modification of the expression, flipping the sign of the term 
containing ^1 , may yield the correct answer. Our numerical value can be used, in a scheme proposed 
by Kotikov, Lipatov and Velizhanin, to estimate the two constants in the strong-coupling expansion 
of the cusp anomalous dimension that are known from string theory. The estimate works to 2.6% 
and 5% accuracy, providing non-trivial evidence in support of the AdS/CFT correspondence. We 
also use the known constants in the strong-coupling expansion as additional input to provide 
approximations to the cusp anomalous dimension which should be accurate to under one percent 
for all values of the coupling. When the evaluations of the integrals are completed through the 
finite terms, it will be possible to test the iterative, exponentiated structure of the finite terms in 
the four-loop four-point amplitude, which was uncovered earlier at two and three loops. 

PACS numbers: 11.15.Bt, 11.15.Pg, 11.25.Db, 11.25.Tq, 12.60.Jv 
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I. INTRODUCTION 



Maximally supersymmetric A/" = 4 Yang-Mills theory (MSYM) has attracted a great 
deal of theoretical interest over the years. It is widely believed that the 't Hooft (planar) 
limit of MSYM, in which the number of colors is taken to infinity, is dual at strong 
coupling (A = g'^Nc oo) to weakly-coupled gravity in five- dimensional anti-de Sitter 
space [II] . The duality implies that the full quantum anomalous dimensions of various series 
of gauge-invariant composite operators are equal to energies of different gravity modes or 
configurations of strings in anti-de Sitter space , |6|. Heuristically, the Malda- 

cena duality conjecture hints that even quantities unprotected by supersymmetry should 
have perturbative series that can be resummed in closed form. The strong-coupling limits 
of these resummed expressions, possibly supplemented by non-perturbative contributions, 
should match results for the appropriate observables in weakly-coupled supergravity or string 
theory. 

This intuition does appear to apply to the higher-loop on-shell scattering amplitudes of 
color-non-singlet gluons, even though the Maldacena conjecture does not directly address 
on-shell amplitudes of massless colored quanta. There is now significant evidence of a very 
simple structure in the planar limit. In particular, the planar contributions to the two- loop 
and three-loop four-gluon amplitudes have been shown to obey iterative relations p, 18[ : The 
dimensionally regularized amplitudes {d = 4 — 2e) can be expressed in terms of lower-loop 
amplitudes, along with a set of three constants for each order in the loop expansion. 

An analogous relation is conjectured to hold for generic maximally-helicity-violating 
(MHV) n-gluon scattering amplitudes, to all loop orders [3, B]. The MHV amplitudes are 
those for which two of the gluons have negative helicity and the remaining (n — 2) have pos- 
itive helicity, or the parity conjugate case. Indirect evidence for the extension to the n-point 
MHV case was provided first by studying the consistency of collinear limits at two loops j3] ■ 
Later, the iteration relation was demonstrated to hold directly for the two-loop five-gluon 
amplitude, for the "even" terms in the am plit ude ol and soon thereafter for the "odd" terms 
as well, i.e. for the complete amplitude 10|. (Odd and even refer to the behavior of the 
ratio of terms in the loop amplitude to the tree amplitude, under parity.) After infrared 
divergences have been subtracted from the proposed all-loop iterative relation, the resulting 
finite remainder is neatly proportional to an exponential of the product of the one-loop fi- 
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nite remainder with the so-called cusp anomalous dimension. Presumably the weak-strong 
duality between anti-de Sitter space and conformal field theory (AdS / CFT) plays a role in 
this simplicity. 

The form of the proposed iterative structure of multi-loop planar MSYM is based on 



the understanding o 
in gauge theory [ll 



how so: 
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t and col 
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mear m 
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13, 



rared singularities factorize and exponentiate 



18| . For the pole terms in the amplitudes, 



such behavior is universal, and holds in any massless gauge theory. What is remarkable in 
planar MSYM is that the finite terms in the MHV scattering amplitudes can be organized 
into the same exponentiated form as the divergent terms. In a non- super symmetric theory. 



QCD, exponentiation of finite terms has also 



3een observed in the context of threshold 



resummation of the Drell-Yan cross section 



19| . In the case of four-gluon amplitudes in 



MSYM, the behavior is valid for arbitrary values of the scattering angle (ratio of t/s). For 
amplitudes with more than four gluons, there are many kinematical variables, and so the 
constraints imposed by the iterative structure are even stronger. It is clearly of interest to 
test whether this structure persists beyond three loops. In this paper we shall provide an 
integral representation for the planar four-loop four-gluon amplitude. Our result will enable 
such a test to be performed at four loops, once the relevant integrals have been evaluated 
to sufficiently high order in e. 

Another remarkable property of planar MSYM is the integrability of the dilatation opera- 
tor, interpreted as a Hamiltonian, for many sectors of the theory. Integrable structures were 
identified in anomalous dimension matrices in QCD a while ago [20]. In planar MSYM, Mi- 



nahan and Zarembo 



2l[ | mapped the one-loop dilatation operator for non-derivative single- 



trace operators to an integrable spin-chain Hamiltonian, and used a Bethe ansatz to compute 
the anomalous dimensions. Such integrable structures have since been extended to higher 
perturbative orders for various sectors of the theory 



22 



23 



2J]. They are also known to 



be present at strong coupling, from the form of the classical sigma model on target space 
AdSg X S"^ [25|. (Berkovits has given a formal argument that the integrability extends to the 



quantum level on the world sheet [26|].) The iterative structure of MSYM amplitudes may 
somehow be related to integrability. If an infinite number of conserved charges are present, 
the form of the quantum corrections should be severely constrained, as it would be by the 
proposed iterative structure j?, 8|. An iterative structure may also arise in correlation func- 
tions of gauge- invariant composite operators in planar MSYM 27(1 ; but its precise structure. 
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if it exists in this context, has not yet been clarified. 

Integrabihty is a powerful computational tool. Integrability, or in some cases, the as- 
sumption of integrability, has been employed to compute a variety of one-loop and multi-loop 



anomalous dimensions in planar MSYM from Bethe ansatze 



21 



22 



231] and from the Baxter 



equation 28|, 1291]. One of the most interesting developments along these lines has been the 



all-orders proposal of Eden and Staudacher (ES), based on an asymptotic all-loop Bethe 
ansatz [s^ , for the large-spin limit of the anomalous dimensions of leading-twist operators 
in MSYM jsit. This quantity, 7A'(as), is also referred to as the cusp (or sometimes, soft) 
anomalous dimension associated with a Wilson line. Equivalently, it represents the leading 
large-x behavior 32| of the DGLAP kernel Pii{x) for parton evolution, i ^ i, 



+ i?(as) 5(1 — x) + . . . , as X 1. 



2(l-x), 

Taking Mellin moments, 7(j) = — dxx^~^P{x), we see that the cusp anomalous dimension 
gives the dominant behavior of the leading- twist anomalous dimensions as the spin j oo, 

1 



7(j) = 7T7x(«s) Hj) + 0{f). 



1.2 



From the asymptotic all-loop Bethe ansatz, ES derived an integral equation for a fiuctu- 
ation density a, from which the cusp anomalous dimension can be determined. The integral 
equation is straightforward to solve perturbatively in a^. In terms of the expansion param- 
eter 



87r2 



27r 



the ES prediction for (one quarter of) the cusp anomalous dimension is 31 ] 

lK{a) 
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Here fnfa) can be identified with the first of a series of three constants (per loop order) 
appearing in the iterative relation for the four-gluon amplitude P, |8|. The first three terms 



of eq. (11.61) agree with previously-known results |8|, l33|, l3J, l35|, l36| , so the new predictions 
begin with the a'* term. The conjecture (11.51) has also been arrived at recently by Belitsky, 
using a proposed generalization of the Baxter equation to all loop orders y^^i- 

In QCD the three-loop cusp (soft) anomalous dimension has been computed by Moch, 
Vermaseren and Vogt as part of the impressive computation of the full leading-twist anoma- 



lous dimensions 



37| needed for next-to-next-to-leading order evolution of parton distribu- 



tion functions. (Terms proportional to the number of quark flavors were obtained first in 

n 

ref. |38|.) Kotikov, Lipatov, Onishchenko and Velizhanin (KLOV) [35] made the inspired 
observation, based on evidence at two loops [39], that the MSYM anomalous dimensions 
may be obtained simply from the "leading-transcendentality" contributions in QCD. The 
cusp anomalous dimensions are polynomials in the Riemann ( values, (n = C{n), or their 
multi-index generalizations, Crai.n2,...- (These cannot show up below five loops.) The degree 
of transcendentality of is just n, and the transcendentality is additive for products of 
( values. At L loops, the leading-transcendentality contributions to the cusp anomalous 
dimension have degree (or weight) equal to 2L — 2. All MSYM leading-twist anomalous 
dimensions computed to date have had uniform, leading transcendentality, whereas the cor- 
responding QCD results contain an array of terms of lower transcendentality, all the way 
down to rational numbers. 

The KLOV conversion principle applies to the leading-twist anomalous dimensions for 
any spin j, with an appropriate definition of leading transcendentality for the harmonic 
sums Srfi{j) that appear in the results. Using assumptions of integrability, Staudacher 



confirmed the three- loop KLOV result through j = 70 40|, |4l|], building on earlier work 



of Beisert, Kristjansen and Staudacher [22| at j = 4. Eden and Staudacher extended this 
analysis to the three- loop cusp anomalous dimension (the j —>■ oo limit), in the course of 
arriving at their all-orders proposal (I1.5P based on integrability The three- loop cusp 
anomalous dimension in MSYM was independently determined from the 1/e^ pole in the 
three-loop four-gluon scattering amplitude , providing a confirmation of the KLOV result 
in the limit j — oo. An important feature of the ES proposal ( 11. 5p is that it is consistent 
with KLOV's observation that the MSYM anomalous dimensions are homogeneous in the 
transcendentality, at least through three loops. 

The ES proposal emerges from mapping single-trace gauge-invariant operators to spin 
chains. The dilatation operator, whose eigenvalues are anomalous dimensions, is mapped 
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to the spin-chain Hamiltonian. The form of the S matrix for this spin chain is fixed, up 
to an overall phase, called the dressing factor {3], by the superconformal [PSU(2, 2|4)] 
symmetry of both AdSs x and Af = 4 supersymmetric gauge theory. At one and two 
loops, superconformal symmetry in combination with explicit calculations fixes the dressing 
factor to be 1. At higher loops, it is not so constrained. A nontrivial dressing factor is 



required by the strong-coupling behavior 



42 



43 



44| . and there have been several recent 



45 



46|. Yet the 



investigations of it using properties such as worldsheet crossing symmetry 
order in the weak-coupling expansion at which it becomes nontrivial is still uncertain. 

The presence of a dressing factor at three loops would lead to a shift in the anomalous 
dimension at four loops. For example, ES have proposed 3l|] a modification of the asymptotic 
Bethe ansatz SOj at this order, with coefficient /3, which is consistent with the presently 
known integrable structure. This modification alters the predicted four-loop anomalous 
dimension in eq. f 1 1.6 1) to 

tt'^-CI + s/^Cs (1.7) 



73 
2520 

Thus a calculation of the four-loop cusp anomalous dimension has the potential to probe a 
nontrivial dressing factor. In a new preprint 47|, Beisert, Eden and Staudacher (BES) have 
shown how to extend the above dressing factor to all orders in the coupling, as well as to 
ensure its consistency with other constraints. This leads to the prediction that /5 = Cs- 

In this paper, we perform this calculation, in order to test whether the ES all-orders 
proposal gives the correct result, or to reveal how it needs to be modified if it doesn't. We do 
so by evaluating the infrared poles of the planar four-loop four-gluon scattering amplitude 
through the order at which the four-loop cusp anomalous dimension appears. The 

form of the infrared singularities at all loop orders are fully understood [15,], up to a set of 
constants that must be computed explicitly. At this undetermined constant is precisely 
the cusp anomalous dimension appearing in the ES formula. 

First, though, we need a representation of the planar four- loop four-gluon amphtude. 



Rather than construct t 
unitarity method 



48 



lis representation from Feynman diagrams, we shall employ the 



49 



50 



51 



two- and three-loop amplitudes 



52 1 . This method was also used to construct the planar 
53| . Because the unitarity method builds amplitudes 



at any loop order from on-shell lower-loop amplitudes, structure uncovered at the tree and 
one-loop levels can easily be fed back into the construction of the higher-loop amplitudes. 
We will find that the planar four-loop amplitude can be expressed as a linear combination 



of just eight four-loop integrals. 

A very interesting feature of the integrals appearing in the planar four-point amplitudes 
through four loops is that they are all, in a well-defined sense, conformally invariant. To 
analyze the conformal properties of potential four-loop integrals, we make use of the recent 
description of such integrals by Drummond, Henn, Sokatchev, and one of the authors 54l |. 

Once we know what four-loop four-point integrals enter into the amplitude, we must 
compute them explicitly through 0(e~^). Here we make use of important recent advances 



61 



62l |. In particular, we use the program 



in multi-loop integration [55|, 156|, 157|, 158|, 159|, 160|, 

MB [g^ to automatically extract poles in e from the Mellin-Barnes representation of loop 
integrals, and to integrate the resulting contour integrals. We carry out the integrations 



analytically for coefficients of the first five poles, through 1/e'^. These coefficients 



are expressed in terms of well-studied functions, harmonic polylogarithms (HPLs) [63 



making it straightforward to confirm the expected infrared structure for arbitrary values of 
the scattering angle. For the coefficients of the Xje" and poles, our analysis is numerical. 
We evaluate the amplitude at four kinematic points, (s,t) = (—1, —1), (—1, —2), (—1, —3), 
and (—1,-15). Numerical evaluation suffices because the expected behavior is completely 
specified at order 1/e^, and specified up to one unknown, but predicted, constant, /q^"*, at 
order 1/e^. We find consistent results from all four kinematic points. 

We also need to evaluate the infrared singular terms to the same order 1/e^. These may 
be expressed in terms of lower-loop amplitudes. For this purpose, we must expand the one-, 
two-, and three-loop amplitudes to 0{e^^, C^(e^) and 0(e°), respectively. Fortunately, these 
are precisely the orders required to test the full iteration relation at three loops, so the needed 
lower- loop expressions can all be found, in analytic form, in ref . Q| . (These results are based 
partly on the earlier evaluation of the double-box Issl] and triple- ladder 57 1 integrals.) 

We find that the ES conjecture is incorrect, although our numerical results suggest that 
a simple modification of the four-loop ES prediction may yield the correct answer. In 
particular, if we fiip the sign of the term in their prediction, we obtain a value within the 
error bars of our result. If we choose to interpret the modification as a dressing factor within 
the form taken in ref. [3l|], it suggests taking the value /3 = for their parameter. This 
value would violate the apparent "uniform transcendentality" observed to date generally for 
quantities in the A/" = 4 theory, for example in the anomalous dimension of the higher-twist 
operator Tr(X^Z^) [311 • The same value, /? = Cs; was suggested independently by 
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BES [47] , based on properties of the spin-chain model. This coincidence, modulo caveats we 
shall discuss in section IVIIt provides some support to violation of uniform transcendentality 
in quantities other than the cusp anomalous dimension. This possibility can be tested via 
other perturbative computations; if uniform transcendentality is nonetheless maintained, 
our result mi ght instead imply that the form postulated for the dressing factor in ES and 



elsewhere 



42 



43 



65, 



66 



671 . l68j is not general enough. 
We also use our four-loop results to investigate the strong-coupling limit of the cusp 
anomalous dimension. In the AdS/CFT correspondence, this quantity can be computed 



from the energy of a long, folded string, spinning in AdS^ 



69| . The first two coefficients in 



the strong-coupling, large- A'c limit of the cusp anomalous dimension 



lave been determined 



using a semi-classical expansion based on this string configuration [33 



69 



70 



7l|. We shall 



discuss how our four-loop result can be used to give a remarkably accurate estimate for 
these coefficients. We employ an approximate formula devised by Kotikov, Lipatov and 
Velizhanin 



34 



35j to interpolate between the weak- and strong-coupling regimes. Using 



our four-loop result as input to this formula, the first two strong-coupling coefficients are 
estimated to within 2.6% and 5%, respectively, of the values computed from string theory. 
This agreement provides direct evidence in support of the AdS / CFT correspondence as well 
as a smooth transition between weak and strong coupling. 

An even better approximation for the cusp anomalous dimension can be found by incor- 
porating into the interpolating formula the string predictions for the first two coefficients 
in the strong-coupling expansion. Based on our success in accurately estimating the two 
leading strong-coupling coefficients, we expect this improved approximation to be accurate 
to within a few percent, for all values of the coupling. Curiously, our approximate formula 
predicts that the third coefficient in the strong-coupling expansion should be very small, and 
may even vanish. The formula also predicts the numerical value of the five-loop coefficient. 
This value turns out to be extremely close to a simple modification of the five-loop ES pre- 
diction, fiipping the signs of the terms containing odd ( values, as at four loops. We have 
confirmed our analysis using Fade approximants, which also give insight into the complex 
analytic structure of /o(a). 

Our representation of the planar four-loop four-gluon amplitude in terms of eight four- 
loop integrals can be used for more than just the extraction of the cusp anomalous dimension. 
As mentioned above, it can also be used to check the proposed iterative structure at four 
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loops. In order to do so, one would need to evaluate the integrals all the way through the 
finite terms, C(e°), instead of just the level carried out in this paper, C(e~^). One would also 
need to evaluate all integrals appearing in the lower-loop amplitudes to two orders higher 
in e. 

This paper is organized as follows. In section [Tll we review the iterative structure of 
MSYM loop amplitudes, commenting in particular on how the cusp anomalous dimension 
appears in the infrared singular terms. In section IIIIl we present the construction of the 
four-loop amplitudes via the unitarity method and also describe the conformal properties of 
the resulting integrals. In section IIVI we give analytical results for the amplitudes through 
0{e~^) and numerical results though 0(e~^), allowing us to extract a numerical value for the 
four-loop cusp anomalous dimension. The four-loop anomalous dimension is then used in 
section IVIII to estimate the coefficients that appear at strong coupling and also to estimate 
higher-loop contributions to the cusp anomalous dimension. Our conclusions are given in 
section rVlIII Two appendices are included, one presenting Mellin-Barnes representations for 
the integrals appearing in the four-loop amplitude, and one reviewing properties of harmonic 
polylogarithms. 



II. ITERATIVE STRUCTURE OF MSYM LOOP AMPLITUDES 



In this paper we consider the planar contributions to gluonic scattering in MSYM with 
gauge group SU{Nc), that is, the leading terms as Nc oo. We do not discuss subleading- 
color contributions; at present they do not appear to have a simple iterative structure 

The leading-color terms have the same color structure as the corresponding tree ampli- 
tudes. The leading- Ai'c contributions to the L-loop SU{Nc) gauge-theory n-point amplitudes 
may be written as, 

L 

J2 Tr(TMi)...TM.))AW(a(l),a(2),...,aH), (2.1) 

where 7 is Euler's constant, and the sum runs over non-cyclic permutations of the external 
legs. In this expression we have suppressed the (all-outgoing) momenta fcj and helicities Aj, 
leaving only the index i as a label. This decomposition holds for all particles in the gauge 
sup er- mult iplet, as they are all in the adjoint representation. The color-ordered partial 
amplitudes An are independent of the color factors, and depend only on the kinematics. For 
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,n-2 



(47r)2-^ 



MSYM, supersymmetric Ward identities 72j imply that the four-gluon hehcity amphtudes 
(++++) and ( — (as well as their parity conjugates) vanish identically. Furthermore, 
the nonvanishing four-point (MHV) amplitudes are all related by simple overall factors. 
Hence we do not need to specify the helicity configuration, i.e. whether the color ordering 

is ( ++) or (-+-+). 

It is convenient to scale out a factor of the tree amplitude, and work with the quantities 
Mi^^ defined by 

Mi'\p;e)=Al!^\p)/A^^\p). (2.2) 

Here p indicates the dependence on the external momenta, p = {si2, S23, ■ ■ .}, where = 
(fcj+fcj-Li )^ are invariants built from color-adjacent momenta. The iteration relation proposed 



in ref. 



then takes the form. 



L=l 



Mnip)^l + J2 ^'Mn\p: = exp J2 (f\e)M!,'\p; le) + C« + E«(p; 6 



1=1 



(2.3) 



In this expression, the factor, 

a^^(4vre--)% (2.4) 

keeps track of the loop order of perturbation theory, and coincides with the prefactor in 
brackets in eq. (12. ip . (It becomes equal to a in eq. ffL3|) as e 0.) The quantity Ml^\p] le) 
is the one-loop amplitude, with the tree amplitude scaled out according to eq. (12. 2p . and 
with the substitution e ^ le performed. (That is, it is evaluated in d = 4 — 2/e.) Each f^^\e) 
is given by a three-term series in e, beginning at 0{e^), 

f'\e) = f^^ + ef?+e'fl\ (2.5) 

The objects ff^, A; = 0, 1, 2, and C(') are pure constants, independent of the external kine- 
matics p, and also independent of the number of legs n. We expect them to be polynomials 
in the Riemann values Cm with rational coefficients, and a uniform degree of transcenden- 
tality, which is equal to 2/ — 2 + A; for and 21 for C'^^\ (Multiple zeta values Cni,n2,... rnay 
also appear, but there are no independent ones of weight less than eight, and so they can 
only appear starting at five- loop order.) The /^'■' and C*^'^ can be determined by matching 
to explicit computations. The one-loop values are defined to be, 

/«(e) = l, C«=0, £;W(p;e)=0. (2.6) 
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The E'i\p-e) are non-iterating 0{e) contributions to the /-loop amphtudes, which vanish 
as e — * 0, En\p;0) = 0. These terms contribute to the exponentiated form of the am- 
phtudes (12.31) even for e ^ because they can appear multiphed by the infrared-divergent 
parts of the one-loop amplitude Mn\p', It). After canceling the infrared divergences between 
real emission and virtual contributions, such terms should not contribute to infrared-safe 
observables. 

The first two values of /^'^(e) in the three-term expansion (1231) . namely f^'^ and ff\ can 
be identified with quantities appearing in the resummed Sudakov form factor 

1 



Ai) 
Jo 



n 



(0 



(0 

Tk : 

Ml) 



The first object, /q''', is identified with the /-loop cusp anoma' 



/q'"* and Qq' = (2/1) ff are known through three loops |8l. |33 



(2.7) 
(2.8) 

ous dimension. The quantities 



35 



36 



39|, 



Jo 

f{2) 

Jo 

A3) 

Jo 



-C2, 



(2.9) 



(see also eq. (11.61) ). and 



Ml) 

yo 

M2) 

yo 



yo 



(3) 



-C3, 

4C5 + y C2C3 • 

c(4) 



(2.10) 



A principal task of this paper is to compute /q and compare the result with the predic- 
tion (ITHD . 



„ n 

Equation (12.31) is equivalent [8|] to 

MI^Hp; e) = X(^)[M«(p; e)] + f'^\e) M'^\p-, Le) + C^^) + Ei^\p; e) 



where the quantities Xn^'^ = X^^'[Mn'] only depend on the lower-loop amplitudes Mn'{p; e) 
with / < L. The xi^^ can be computed simply by performing the following Taylor expansion, 



<L)nAi) 



(2.11) 



(0/ 



Xf) [Mf] = M(-)- In 1 + J] a'M« 



1=1 



(2.12) 
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Eqs. f l2.1ip and fl2.12p express the L-loop amplitude explicitly in terms of lower-loop ampli- 
tudes, plus constant remainders. Here we need the values of Xn^^ for L = 2, 3, 4, 

,1 2 



X(3)[M«] 



MP 



1 3 



MP 



1 4 



'Mi^) + M^mP + 



(2.13) 
(2.14) 
(2.15) 



We note that the exponentiated result (12.31) leads to a simple exponentiated form for 
=u.tably-deaned "finite .en,ainde..' ' a^sooated w.th the n,uM-loop a,npMude= We 

define 

Fi.'\p: e) = M(^) - J2 ^^"'^ ^1'^ ' (2-16) 



1=0 



where the ip^''\p; e) are iteratively-defined divergent terms, and mP = 1. After some alge- 
bra, one finds that Ip\p] e) and Fn^\p; e) obey iterative relations very similar to eq. (12.111) . 
In the limit as e — 0, the relation for F^\p] e) becomes 



(2.17) 



Because e has disappeared from eq. (12.171) . it can be solved neatly for Fn^\p; 0) for any L, 
in terms of the o„e-loop remainder F»'(p; 0) alorre 1^. The solution can be represented as, 



OO p OO 

^„(p;0) = l + ^a^i^f)(p;0) = exp a' (/o^M'Hp; 0) + C« 



L=l 



exp 



1=1 



-j^(a) FP{p;0) + C{a) 



(2.18) 



where C{a) = Xli^i ^''-'-'a' and we used the relation (12. 7p of Jq^ to the cusp anomalous 
dimension. The result for F'P^ (p; 0) is given by the term in the Taylor expansion of the 
exponential. 

Next we present the specific forms of the iterative amplitude relations (12. lip through four 
loops, specializing to = 4. The two-loop version is |3] 

Mf (p; e) = ^ [Mi'^ (p; e)]\ f^'^ (e) M^'^ (p; 2e) + C^'^ + 0(e) , (2.19) 



where 



(2.20) 
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and the constant C*^^-* is given by 



2^2- 



(2.21) 



The three-loop version, exphcitly verified in ref. 

3 



IS 



Mf (p;e) 



Mi'^ (p; e) J V Mi^) (p; e) Mf (p; e) + (e) Mi^^ (p; 3 e) 



where 



/ 



(3)/ 
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C4 + e(6C5 + 5C2C3) + e'(ciC6 + C2CI) 



and the constant C^^-* is given by 



341 2 \ / 17 2 \ 2 
^ +gCljC6+(^-y + gC2jC3• 



(2.22) 



(2.23) 



(2.24) 



The constants Ci and C2 are expected to be rational numbers. They drop out from the right- 
hand side of eq. fl2.22p because of a cancellation between /g^"* and C^^K A computation of 
the three-loop five-point amplitude, or of the three-loop splitting amplitude, could be used 
to determine them. 

The four-loop iteration relation would have the following form, 
I r 1 4 r 1 2 



(p;6) 



4 



+ 



Mr(p;e 

1 
2 



m1^) (p; e) M^' (p; e) + ^ (p; e) Mf (p; e 



r(2) 



Mf (p; 6) ' + /W(e) Ml^'(p; 4e) + C^^) + 0{e) 



(1), 



(2.25) 



As we shall not be computing the 1/e and finite terms in the present paper, we cannot do 
more here than verify the (universal) divergent terms and extract the value of the four- 
loop cusp anomalous dimension. We leave to future work the important task of verifying 
eq. (12.251) . using the integral form of the four-loop amplitude presented in this paper. 



III. CONSTRUCTION OF FOUR-LOOP PLANAR MSYM LOOP AMPLITUDE 



The unitarity method 



48 



49 



5G, 



51 



52l | is an efficient way to determine the represen- 



tations of loop amplitudes in terms of basic loop integrals. The coefficients of the loop 
integrals are obtained by sewing sets of on-shell tree amplitudes. If we are using a four- 
dimensional form of the unitarity method, the tree amplitudes can be significantly simplified 
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before sewing. At one-loop supersymmetric amplitudes are fully determined from their four- 
dimensional cuts, but unfortunately for higher loops no such theorem has been proven. In 
the present calculation, we will therefore use D-dimensional unitarity, for which we will need 
the tree amplitudes to be evaluated without assuming the four-dimensional helicity states 
for the external legs. These amplitudes are nonetheless simpler than the completely off-shell 
amplitudes that would implicitly arise in a conventional Feynman-diagram calculation. For 
MSYM, a key feature is that the on-shell tree amplitudes have the full Af = 4 supersymmetry 
manifest, in the form of simple S'-matrix Ward identities [72]. It is impossible to maintain 



the full J\f = 4 supersymmetry in any off-shell formalism, 
imply the equations of motion via the Bianchi identities 



3ecause the superspace constraints 



3- 



The unitarity method derives 



its efficiency from the ability to use simplified forms of tree amplitudes to produce simplified 



loop integrandS; (To maintain the supersymmetry, we apply the 



bur-dimensional helicity 



(FDH) scheme 7^ , a variation on dimensional reduction (DR) 75| , in performing the sum 
over intermediate gluon polarization states.) 

The unitarity method expresses the amplitude in terms of a set of loop integrals. In 
general gauge theories, such as QCD, the number of required integrals proliferates rapidly 
as the number of loops increases, and sophisticated algorithms based on integration-by- 



parts identities [76l have been devised to relate such integrals to a smaller class of master 
integrals, successfully through two loops |77| . Fortunately, the number of required integrals 
grows much more slowly for gluon-gluon scattering in planar A/" = 4 super Yang-Mills theory. 
At L = 1, 2, 3 the respective numbers are 1, 1, 2, and the required integrals are all shown in 
fig. [H We will see that at L = 4 eight integrals are required. 



The result for the one-loop four-point amplitude is [781] 



(3.1) 



2 

where the Mandelstam variables are s = {ki + and t = {k2 + /cs)^. The factor of 1/2 
in eq. (13.11) follows from our normalization convention for Ait\ which is defined by eq. (12.11) . 
The one- loop scalar box integral I^^\s, t) (multiplied by a convenient normalization factor), 
depicted in fig. [H is 

I^'\s,t) = stli'\s,t), (3.2) 

where I^^ (s, t) is defined in eq. (Bl) of ref. . We absorb the factor of st into the definition of 
the integrals we use, because it cancels a factor of 1/ (st) appearing in the explicit expression 
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1 4 
(1) 



(2) 



(3)a 



(3)b 



FIG. 1: Integrals required for gg gg scattering in planar MSYM at one loop (1), two loops (2) 
and three loops ((3)a and (3)b). The box (1), planar double box (2) and three-loop ladder (3)a 
integrals are scalar integrals, with no loop-momentum dependent factors in the numerator. The 
tennis-court integral (3)b contains a factor of {li + 12)'^, where li and I2 are marked with arrows in 
the figure. 

for I^\s,t), and matches the form in which it appears in the MSYM amplitudes. This 
integral is given in terms of HPLs in eq. (Bl) of that reference, through the order we 
require here, O(e^). The higher-order terms in e are needed to be able to evaluate terms in 
the infrared/iterative representation (12.251) of the four-loop amplitude through (9(e^^). For 
example, in the term [M4^'*(e)]^ oc [X''"'^^(e)]^ in eq. (12.251) . if one takes the leading term 
from three of the four factors, then the coefficient of in the fourth factor contributes to 
the (9(e~^) term in the product. 



The planar two-loop MSYM four-point amplitude is given by 53 1 



Mf (e) 



The two-loop scalar double-box integral, shown in fig. [T], is 



(3.3) 



(3.4) 



where J4 (s, t) is defined in eq. (B4) of ref. ^]. As at one loop, we have rescaled the integral 
to remove the rational prefactor. This integral was first evaluated through 0(e°) in terms of 



polylogarithms [55|. Because the infrared/iterative expression (I2.25P contains, for example, 
[MP(e)]2 oc [J(2)(e)]2, and because the expansion of I^'^\e) begins at order we need 
its expansion through 0{e^). This expansion is presented in terms of HPLs in eq. (B5) of 
ref. Is]. 

The three-loop planar amplitude is given by 

1 



53j 



Mf (e) 



J(3)"(s, t) + 2 J(=^)^(t, s) + l^^^\t, s) + 2 J(=')^(s, t) 



(3.5) 
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The scalar triple-ladder and non-scalar "tennis-court" integrals, illustrated in fig. [T], are 

J(=^)" = st2/P(s,t), (3.6) 



where I^^^{s,t) and I^^^{s,t) are defined in eqs. (3.1) and (3.2), respectively, of ref. 
Because these integrals multiply X*^^-' in the term M^\e)M^\e) in eq. fl2.25p . we need their 
expansion through 0{e^). These expansions were first carried out in terms of HPLs for I^^^^ 



in ref. 57|], and for I^^^ in ref. 8(]. The results are collected in eqs. (B7) and (B9) of ref. 



The coefficients of the integrals in the two- and three-loop expressions 03.31) and (13. 5p 



were originally determined [53|] using iterated two-particle cuts. Such cuts can be evaluated 
to all orders in e because Af = 4 supersymmetry relates all nonvanishing four-point ampli- 
tudes; therefore precisely the same algebra enters as at one loop, for which it leads to the 
amplitude (13. ip . More generally, an ansatz for the planar contributions to the integrands 
was proposed in terms of a "rung insertion rule" [3 S] to be described below, which was 
based largely on the structure of the iterated two-particle cuts. At three loops, the planar 
integrals generated by the rung rule can all be constructed using iterated two-particle cuts. 
Also, the three-loopplanar amplitude (13. 5p has the correct infrared poles and a remarkable 
iterative structure [8|, so there is little doubt that it is the complete answer. 

However, beyond three loops — and even at three loops for non-planar contributions — 
the rung rule generates graphs that cannot be obtained using iterated two-particle cuts. It 
is less certain that the rung rule gives the correct results for such contributions. Indeed, 
we shall see that there are additional, non-rung-rule, contributions to the planar amplitude 
beginning at four loops. 

Nevertheless, we start constructing the planar four-loop MSYM amplitude using the 
diagrams generated by the rung rule. According to this rule, each diagram in the planar 
L-loop amplitude can be used to generate planar (L-|- l)-loop diagrams as follows: First, one 
generates a set of diagrams by inserting a new line joining each possible pair of internal lines. 
Next, one removes from this set all diagrams with triangle or bubble subdiagrams. Besides 
the scalar propagator associated with the new line, one also includes an additional numerator 
factor for the diagram, beyond that inherited from the L-loop diagram, of + /2)^- Here h 
and I2 are the momenta fiowing through each of the legs to which the new line is joined. Each 
distinct (L -|- l)-loop contribution is counted once, even if it can be generated in multiple 
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FIG. 2: "Rung-rule" contributions to the leading-color four-loop amplitude, in terms of integral 
functions given in eqs. ()A9p - ()A14p . An overall factor of st has been suppressed in each figure, 
compared with the definitions in eqs. ()A9p - (jA14p . 

ways. (Contributions corresponding to identical graphs but with different numerator factors 
should be counted as distinct.) Rung- rule diagrams have also been referred to as "Mondrian 
diagrams" because of their visual similarity j8|. 

At one loop, the only no-triangle graph for the four-point process is the box graph depicted 
in fig. [U In going to two loops, there is only one inequivalent way to add a rung to the 
box graph without creating a triangle. Connecting opposing sides of the box, we form the 
planar double box in fig. [TJ (One might also imagine attaching a propagator between two 
adjacent external legs. However, this operation yields the same double box.) To generate 
the three-loop graphs, we add a rung, either vertically or horizontally, to one of the boxes 
of the two- loop box diagram, thus yielding the triple ladder (3) a and tennis-court integral 
(3)b of fig. dJ (Attaching a propagator between external legs again gives nothing new.) Of 
course, there are several permutations of external legs present for any given type of graph; 
all of them are produced by the rung rule. Here we are identifying only the topologically 
distinct integrals that arise. 

What happens when we try to add rungs to the three-loop integrals? There are two 
inequivalent ways to add a rung inside either the left- or the right-most box of the triple 
ladder integral in fig. [U these give the integrals of fig. [2](a) and (c). Adding a vertical rung 
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(ii) 



(iii) 




il ImI til 



(V) 



(vi) 



FIG. 3: Generalized cuts that provide information about the planar four-loop amplitude, (i) A 
two-particle cut separating a tree amplitude from a three-loop amplitude, (ii) A two-particle cut 
separating a one-loop amplitude from a two-loop amplitude, (iii) A "3-3" cut separating the 
amplitude into a product of three tree amplitudes, (iv) An "upper-2-3-lower-2" cut separating 
the amplitude into a product of four tree amplitudes, (v) A "lower-2-3-lower-2" cut. (vi) A 
"3-lower-3" cut. 

inside the middle box does not yield a topologically distinct integral. Adding a horizontal 
rung inside the middle box yields the integral of fig. [2](d). Inserting a vertical rung inside 
the upper-right box of the tennis-court integral in fig. [T] gives us fig. [2](e), and a horizontal 
one, fig. MJo). Finally, adding a horizontal rung inside the left-side box of the tennis-court 
integral gives us a new kind of integral, shown in fig. Et^f), which has no two-particle cuts. 

The propagators and momentum numerators present in integrals fig. [2](a)-(e) are deter- 
mined by the two-particle cuts, as depicted in fig. [3](i) and (ii). That is, the rung rule is 
guaranteed to be correct for them. There is no such guarantee for integral fig. EJ^f), which 
has no two-particle cut with real external momenta. In order to check it, we need to com- 
pute a three-particle cut. We chose to perform the generalized-unitarity cut of fig. [3](iv), a 
threefold-cut with a central three-particle cut and two secondary two-particle cuts. This cut 
reveals that the rung rule is more robust than might have been expected, based on its origin 
in iterated two-particle cuts: The rule does in fact give the correct form for the numerator 
of the integral in fig. Mj), even though no two-particle cut can detect it. 
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FIG. 4: Non-rung-rule contributions to the leading-color four- loop amplitude, in terms of integral 
functions defined in eqs. (jAlSp and (|A16p . Integral (d2) follows the labeling of integral (d) in fig. El 
and integral (f2) follows the labeling of integral (f). An overall factor of st has been suppressed in 
each figure, compared with eqs. (|A15P and ()A16p . 

Because we have no proof that the (— 2e)-dimensional parts of loop momenta are unimpor- 
tant to the computation, we perform these calculations in D dimensions. For this purpose, 
we need tree amplitudes with (some of) the external states kept in general dimension. In 
the three-particle cuts, one has contributions from three-gluon states, and from states with 
gluons and fermion (gluino) pairs crossing the cut, in addition to states with scalar pairs or 
fermion pairs and a lone scalar. In principle, one could evaluate these cuts by computing all 
the required amplitudes, and summing over the particle multiplet. However, it is easier to 
use a trick. The trick takes advantage of the fact that the MSYM multiplet, consisting of 
the gluon, four Majorana fermions, and six real scalars, can also be understood as a single 
A/" = 1 multiplet in ten dimensions. Instead of summing over the multiplet seen as the 
A/" = 4 multiplet in four dimensions, we sum over the Af = 1 multiplet in ten dimensions. 
This trick reduces the number of types of intermediate states one has to consider, though 
of course the total number of states is unchanged. The loop momenta are in any event kept 
in D dimensions. (The trick is compatible with the FDH regularization scheme 7J|, where 
the momenta are taken to be in D dimensions, but the number of states in the loops is kept 
at the four-dimensional value.) 

The cut of fig. [3](iv) also reveals the presence of a non-rung-rule integral, shown in 
fig. Il](d2). The integral has, of course, no two-particle cuts, but it can be obtained from 
the integral topology of fig. [2](d) by canceling two propagators, labeled 10 and 13. The 
integral of fig. [2](f) could also have been checked using another generalized- unit arity cut, 
the two-fold three-particle cut of fig. [3]^iii). This cut also reveals the presence of the "four- 
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square" integral, shown in fig. Ill^f2). This integral can be obtained from the topology of 
fig. [2t^f) by canceling the propagator labeled 13. The cut of fig. |3t^iii) also detects the inte- 
gral of fig. IKd2). 

Under mild assumptions, which we discuss in section |Vl it is sufficient to compute two 
additional multiple cuts, shown in fig. [3](v) and (vi), in order to rule out any additional 
contributions. We have computed these cuts, and we indeed find that no additional integrals 
appear. 



IV. INTEGRAL REPRESENTATION OF THE FOUR-LOOP PLANAR AMPLI- 
TUDE 

We find that the four-loop planar amplitude is given by 



M, 



(4), 



4 



1 
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J(") (s, t) + J(") (t, s) + 2 JC^) (s, t)+2 J(^) (t, s) + 2 J(^) (s, t) + 2 J(=) (t, s) 
+ (s, t) + J(^) (t, s) + 4 J(^) (s, t) + 4 T^-^) (t, s) + 2 J(^) (s, t) + 2 1^'^ (t, s) 
-2X^'^^\s,t)-21^'^^\t,s)-l'^^^\s,t) \ . (4.1) 

The rung-rule integrals X^'^^ through X*^^) are depicted in fig. [2l The two additional integrals, 
X(d2) Qj^f^ T^^^\ depicted in fig. HI do not follow from the rung rule, and were detected using 
generalized cuts with at least one three-particle channel, as discussed in section IIIII 

The integrals appearing in the four-point amplitude are defined generically, for diagram 

"(x)" by 

/of \f{^) 
d''pd''qd''ud''v (4.2) 
UjPj 

where p, q, u, v are the four independent loop integration variables, and c? = 4 — 2e. The 
product in the denominator of eq. (14.21) runs over the labels of internal lines in the graph (x) . 
(For graphs (b) and (c), line 11 corresponds to a numerator factor, so it should be excluded 
from this product. Similarly, lines 10 and 13 are to be omitted from the denominator product 
for graph (d2), and line 13 from the product for graph (f2).) Each line carries momentum 
Pj, which is some linear function of p,q,u,v and the external momenta. The line label j is 
shown next to each internal line. The numerator factor J\f^^^ is also shown explicitly, to the 
left of the graph for (x). We have omitted an overall factor of st from J\f^^^ associated with 
each integral from the figure, in order to avoid cluttering it. 



21 



For the quadruple-ladder graph (a), and for the non-rung-rule graphs (d2) and (£2), A/"*-^^ 
is completely independent of the loop momentum, and so it may be pulled outside of the 
integral. For the other graphs, at least one factor in J\f^^^ depends on the loop momenta. 
In generating a Mellin-Barnes representation for these integrals, we think of these factors 
as additional "propagators" appearing in the numerator instead of the denominator. Ac- 
cordingly, we attach a line label j to each such factor. The presence of such a factor is also 
indicated graphically by a pair of parallel arrows, marking the lines whose momenta are 
summed, then squared, to generate the numerator factor. 

For example, the quadruple ladder integral (a) is defined by 

J {p — kiY {p — ki — ^2)^ {p — "^"^ — i-. — i-^^^ 



X 



{p — {p — ki — ^2)^ q"^ {p — qY {q — ki — ^2) 

1 

{q — uY {u — ki — k2Y {u ~ {v — ki — ^2)^ {v + 



.e-V-'T / d-pd^'-f f (4.3) 

11^ = 1 Pj 



Similarly, integral (b) is defined by 



X(b) = (_,e-vr"^/^)St^ I d'^pd'^rd\d'v[{v + k,n 

J p^ {p — kiY {p — V — kiY (f- {p — qY {q — V — kiY 

1 



X 



■u^ {q — uY {u — V — kiY {u + k4Yv'^ {v — k2Y {v — k2 — k^Y 
■-.e.^.-^/y/ '^^^ff-f^ff-r (4.4, 



The double arrows indicate that in this case the numerator factor appears squared, {p 

{is+hor = [{v+k,Yr- 



2 \2 

u) 



V. ESTABLISHING THE CORRECTNESS OF THE INTEGRAND 

In this section we justify the result (14. ip for the four-loop planar amplitude, based upon 
our evaluation of the unitarity cuts. Because we have not evaluated all possible unitarity 
cuts, we have to impose some mild assumptions about the types of integrals that should be 
present. We will see that another, stronger assumption of conformal invariance also holds 
for the individual integrals that appear, although we do not require it. In addition, as we 
shall see in the next section, the agreement of the infrared singularities through (9(e^^) with 
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their known form [15|] — up to the one unknown constant at (!?(e^^) — provides a non-trivial 
consistency check on our construction. 

The analysis determining the integrand of the four-loop four-point amplitude proceeds in 
several steps: 

• We assume that there are no integrals with triangle or bubble subgraphs. 

• We classify the four-loop planar integrals of this type topologically. We begin with the 
subset of graphs having only cubic vertices, from which we can obtain the remaining 
graphs. 

• We construct a set of generalized cuts capable of detecting all such integrals. That 
is, each such integral, when restricted to the generalized cut kinematics, should be 
nonvanishing for at least one cut in the set. For it to be nonvanishing, it must have 
a propagator present for each line being cut. We used this set of cuts to deduce the 
terms in the expression (14. ip . Indeed, we find that each such cut of the expression is 
completely consistent with our evaluation of the cut. This step confirms the result, 
under the "no triangle" assumption. 

• Alternatively, we consider the result of assuming that only conformally-invariant in- 
tegrals contribute, when the external legs are taken off-shell so that the integrals be- 
come well-defined (finite) in four dimensions. Such integrals were considered recently 
in ref. {s^. We find that the conformal-invariance assumption is a powerful one; it 
allows all eight integrals contributing to eq. (14.11) to be present, while forbidding all 
but two of the additional no-triangle integrals. The potential contributions of these 
remaining two integrals are easily ruled out by examining the two-particle cuts. 

Next we elaborate on each step in the analysis. 



A. Unitarity construction 

Our assumption, that there are no integrals with triangle or bubble subgraphs in multi- 
loop A/" = 4 super- Yang-Mills theory, is sometimes referred to as the "no triangle hypothesis" . 
(Such an assumption also appears to be applicable to A/" = 8 supergravity, at least at one 
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loop [80].) We now discuss evidence in favor of this hypothesis. First, notice that a bub- 
ble subgraph would lead to an ultraviolet sub divergence. In the absence of cancellations 
between different integral topologies, such sub divergences are forbidden by the finiteness of 
MSYM jsil. Keep in mind that all cancellations between different particles in the supermul- 
tiplet have already been taken into account, so the coefficients of all bubble integrals should 
indeed be zero. 

Next, suppose there were a triangle subgraph in some multi-loop integral. Excise a 
region around the triangle, cutting through open propagators attached to the integral. This 
excised region represents a one-particle-irreducible triangle-type contribution to a one-loop 
n-particle scattering amplitude. If all n cut legs are gluons, then we know that such a 
contribution is forbidden in MSYM, for arbitrary n, by applying loop-momentum power- 



counting to a computation of the one-loop amplitude using background-field gauge 48|, l82 |. 
However, in the present case some of the n legs might be associated with other fields of 
the Af = A supermultiplet. Supersymmetry Ward identities [t^ typically relate many such 
amplitudes to the gluonic case, but we do not know of a general proof. Thus we do not 
claim to have a full proof that all topologies with triangle subgraphs are absent, although 
we strongly suspect that it is the case. 

We now wish to classify the four-loop planar integrals containing no triangle or bubble 
subgraphs. We can perform this classification first for the subset of such graphs that have 
only cubic (three-point) vertices, for the following reason: If a graph contains a quartic or 
higher-point vertex, we can "resolve" the vertex into multiple three-point vertices by moving 
some of the lines attached to the vertex. Such a procedure never decreases the number of 
propagators associated with a given loop. Therefore a no-triangle graph will remain a no- 
triangle graph under this procedure. So, if we know all the no-triangle graphs with only cubic 
vertices, we can get every remaining no-triangle graph by sliding cubic vertices together until 
they coincide, a procedure known as "canceling propagators" . 

The cubic subset can be classified iteratively in the number of loops using the Dyson- 
Schwinger equation. In order to use the Dyson- Schwinger approach for the case at hand, the 
planar four- loop four-point amplitude, we also need to classify the planar cubic no-triangle 
graphs of the following types: one loop for the number of legs n up to 7, two loops for n up 
to 6, and three loops for n up to 5. The result is that there are a total of 13 planar cubic 
four-loop four-point no-triangle graphs. Of the 13 graphs, seven are one-particle-irreducible. 
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FIG. 5: The only one-particle-irreducible purely-cubic four- loop four-point graph with no triangle 
or bubble subgraphs, besides the rung-rule graphs in fig. [2l 

Six of these seven have the topology of the rung-rule graphs shown in fig. [2j The seventh 
is shown in fig. [5l Note that here we are only classifying graphs, and not yet specifying the 
loop-momentum polynomials associated with each graph. 

The six one-particle-reducible graphs can be obtained by sewing external trees onto four- 
loop graphs with either two or three external legs. (There are no planar cubic no-triangle 
graphs at four loops with fewer than two external legs.) From the point of view of the 
cuts, these one-particle-reducible graphs are equivalent to certain of the non-cubic graphs 
obtained by canceling external propagators, so we shall defer their description briefly. 

The next step is to cancel propagators between vertices in the cubic graphs in figs. [2] and 
[5l That is, we merge adjacent three-point vertices into four-point (or higher-point) vertices 
by eliminating the line(s) between them, while insisting on at least four propagators around 
every sub-loop. This procedure is also straightforward to carry out. It generates the two 
non- rung-rule graphs in the expression (14.11) . shown in fig. |U as well as the 16 additional 
graphs shown in fig. [6l We have given each graph a notation which indicates the rung-rule 
graph (or graph (g)) from which it can be generated by canceling one or more propagators. 
For example, graph (ba) is found by canceling propagators 13 and 14 in rung-rule graph (b), 
and graph (ee) is found by canceling propagators 8 and 11 in rung-rule graph (e). 

Some graphs in fig. [6] can be generated from more than one cubic graph. In fact, five of 
the six one-particle- reducible cubic graphs mentioned earlier are equivalent, upon canceling 
external propagators, to diagrams in fig. [6l namely graphs (b2), (da), (ds), (ei), and (gi). 
Hence we do not provide a separate figure for the one-particle-reducible cubic graphs. The 
sixth graph has the form of a massless version of graph (ds), sewn as an external bubble to 



25 



m ^ M ^ m ^ 

(b^) (b2) (b3) (b4) (c^) (di) 
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^ m ^ ^ 

(64) (e^) (e^) (g^) 

FIG. 6: No-triangle planar four-loop graphs, in addition to those given in figs. [2] and [H The 
notation indicates the rung-rule graph in fig. [21 or graph (g) in fig. [5l that a given graph here is 
derived from by canceling propagators. 

a four-point tree graph. Such integrals vanish in dimensional regularization, so we need not 
consider them. 

To complete the direct justification of the four- loop result ( 14. ip . given the no-triangle 
assumption, we just need to show that each graph appearing in figs. Ej HI [5] and [6] can be 
detected by at least one of the generalized cuts we have computed, out of the six cuts shown 
in fig. [3l Table [T] summarizes some of the cuts that detect the no-triangle graphs. For some 
of the graphs, other cuts also detect them, but for brevity they were not listed in the table. 
Almost all graphs appear in multiple cuts. Only two of the graphs, (b4) and (eg), appear in 
a unique cut, the "3-lower-3" cut labeled (vi) in fig. [3l 



B. Conformal properties 

We have finished our justification of the representation (14. ip for the four-loop planar 
amplitude. In the rest of this section we would like to examine the consequences of making 
a stronger assumption than the no-triangle hypothesis. This assumption is that each of 
the integral functions that appears is conformally invariant. Here we are inspired by the 
discussion of conformally-invariant integrals by Drummond, Henn, Sokatchev, and one of 



26 



TABLE I: The no-triangle graphs, and some of the cuts from fig. [3] that detect them. (In some 
cases, the diagram must be rotated or flipped first.) 



Graph 


Cuts 


Graph 


Cuts 


Graph 


Cuts 


(a) 


(i), (ii), (iii) 


(bl) 


(v), (vi) 


(d5) 


(i), (iii), (iv) 


(b) 


(i), (v) 


(b2) 


(i), (vi) 


(ei) 


(i), (iii), (iv) 


(c) 


(i), (ii), (iii), (iv) 


(bs) 


(v), (vi) 


(e2) 


(iii), (iv) 


(d) 


(i), (iii), (v) 


(b4) 


(vi) 


(ea) 


(i), (vi) 


(e) 


(i), (iii), (iv), (vi) 


(ci) 


(i), (iii), (iv) 


(e4) 


(i), (vi) 


(f) 


(iii), (iv), (vi) 


(di) 


(i), (iii), (iv) 


(es) 


(i), (iii), (iv) 


(d2) 


(iii), (iv) 


(da) 


(i), (iii), (iv) 


(ee) 


(vi) 


(f2) 


(iii), (vi) 


(d4) 


(i), (iii), (iv) 


(g) 


(i), (vi) 










(gl) 


(i), (vi) 



the authors |54|. Although the requirement of conformal invariance is natural because of 
the conformal invariance of the theory in four dimensions, we do not have a proof that these 
integrals are the only ones that can appear in the amplitudes. Nevertheless, as we shall see, 
the conformal properties offer a rather useful guide. (It is possible that extensions of the 
conformal-invariance analysis in ref. 831] could be used to prove that only such conformally 
invariant integrals can be present.) 

We actually wish to study the conformal properties of the integrals in four dimensions, 
yet they are ill-defined there because of the infrared divergences associated with on-shell, 
massless external legs. We therefore adopt a different infrared regularization of the integrals 
by taking the external legs off shell, letting kf 0, i = 1,2,3,4, instead of using a dimen- 
sional regulator as in the rest of the paper. We demand that each integral that appears 
be conformally invariant. Actually, the integrals need only transform covariantly, carrying 
conformal weights associated with each of the external legs, in such a way that they can 
be made invariant by multiplying by appropriate overall factors of s or t. In canceling the 
conformal weights using the external invariants, we should not use any factors that vanish as 
the external legs return on shell, in the limit kf —>■ 0. Such factors would lead to power-law 
divergences or to vanishings of the integrals that are too severe, compared to the typical 
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logarithmic dependences on kf from the known form of infrared singularities. Thus integrals 
that require powers of kf to be conformally invariant should not appear in any on-shell 
amplitude. This on-shell restriction turns out to be a powerful one. 

The net result of the conformal-invariance requirement will be that, besides the eight 
integrals already present in eq. (14. ip and figs. [2] and HI remarkably only two other potential 
conformal integrals survive. The first of them is the (ds) graph from fig. O This graph has 
the structure of a potential propagator correction at four loops. However, in the context 
of the gg — > gg scattering amplitude of the Af = 4 theory, it can be excluded very simply, 
without computing any generalized cuts. One only has to use the structure of the three- 
loop amplitude, and the simplest two-particle cut, cut (i) in fig. [3l to see that it cannot be 
present. The second of the additional potential integrals has the topology of integral (d) in 
fig. El but it has a different numerator factor, to be described below. This new integral, (d'), 
is also easy to exclude using the same two-particle cut (i) in fig. [31 

It is rather striking that every integral identified via the unitarity cuts is conformally 
invariant, and that there are only two other conformal integrals, which can be eliminated 
easily via two-particle cuts. Furthermore, the two integrals that are not present differ from 
the eight that are present in how the conformal invariance is achieved, as we shall discuss 
below. 

To analyze the conformal invariance properties, we shall use changes of variables as 



suggested in ref. 5J]. (The same conformal integrals appear in coordinate-space correlators 



of gauge invariant operators 27|]; the coincidence is presumably an accident of there being 
a limited number of conformal integrals.) As an example, consider the two-loop double box 
depicted in fig. [H 



d^p d'^g 



(5.1) 



p^{p - kiY{p -ki- k2Yq^{q - kiY{q - ks - kiY{p + qY ' 

We have taken d = 4, with the ki off shell to serve as an infrared regulator. Next, use the 
change of variables, 

/Cl=X4i, k2=Xi2, ks = X23, k^ = Xu , ^ = ^45, q = Xq^ , (5.2) 

where This choice of variables automatically ensures that momentum is 
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conserved, ki + k2 + + = 0. Note that the external invariants become 

s = {ki + k2f = , t= {k2 + k^f = . (5.3) 
Performing the change of variables ( 15.20 in the double box, we obtain, 

l'^^\xi,X2,X3,X4) = {-ie'^TT~^Yx24xlr^ [ d^Xg d^Xg 2 2 2 2 2 - i^'^) 

I ly-^ ry-^ ly^ ly^ 

J "''45 "''1 5 "''25 -^46 "''36"'' 62"'' 56 

The principal conformal-invariance constraints on integrals constructed from the invariants 
x\j are exposed by performing an inversion on all points, x^jx^. (We cannot impose 

such an inversion on the ki directly, because it would violate the constraint of momentum 
conservation.) Under the inversion, we have 

2 -^ij 14 d X5 4 d X% , , 

^^i^ZIZf' dx5^— ^, dxg^— ^. (5.5) 

X5 Xg 

It is easy to see that the planar double-box integral is invariant under inversion, i.e. 
X*^^^(xi, X2, X3, X4) X^^)(xi, X2, X3, X4). For this result to hold, it is important that the 
unintegrated points xi, X2, X3, X4 appear in the numerator just enough times to cancel their 
appearance in the denominator. The integrated points X5,X6 each appear four times in 
the denominator. The dimensionally-regulated version of the conformally-invariant integral 

X(2) 

(xi, X2, X3, X4) is precisely the form in which the two-loop double box appears in the 
two- loop planar amplitude (13. 3p . 

brmal properties of integrals beyond two loops, it is helpful 



In order to analyze the con 
to follow the discussion of ref. 5^1, and introduce a set of dual diagrams [S^]- We construct 
the dual to a diagrammatic representation of a planar loop-momentum integral by placing 
vertices corresponding to the Xj at the centers of the loops and in between pairs of external 
lines. Denominator factors of xl^ are denoted by drawing dark solid (blue) lines between 
the corresponding vertices. Numerator factors are denoted by drawing dotted lines between 
the corresponding vertices. One solid line crosses each propagator in a loop. The conformal 
weight in each Xj variable is then given by the number of solid lines entering the corresponding 
vertex, less the number of dotted lines. A conformally-invariant integral will have weight 
four at each internal vertex (to balance the weight of the integration measure), and weight 
zero at each external vertex. In the diagrams, we shall omit one dotted line connecting 
external vertices X2 and X4, and another one connecting x\ and X3, in order to simplify the 
presentation. These two omitted lines correspond to the overall factor of st = x\^x\^ omitted 
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FIG. 7: The two-loop planar double box and its dual diagram. The double box is represented 
by light colored lines and the dual diagram by dark (blue) lines. A dark line connecting Xi with 
Xj represents the factor l/xfj. A dotted line signifies a numerator factor of xfj. The momentum 
corresponding to any Xij is given by the sum of momenta of the light lines crossing the dark line 
joining Xj and Xj. An overall factor of st has been removed for clarity. 

from the momentum-space diagrams in figs. [2] and HI Expressions for integrals in terms of 
the Xi variables can be read off quickly from the dual diagrams (and vice versa). 

For example, fig. [7] contains the diagram dual to the double box. Each of the solid 
lines starting at an Xj and ending at an xj corresponds to a factor of appearing in 

eq. (15.41) . while the dotted line corresponds to a factor of = s. With this identification 
the dual figure is in direct correspondence with eq. (15. 4p . after removing one overall factor 
of st = 0:24X^3 (in order to reduce the visual clutter in the diagram). Because the number 
of solid lines minus the number of dotted lines at each of the two internal vertices X5 and 
Xq is four, the integral is conformally invariant with respect to these points. Similarly, since 
each of the external points Xi,X2, x^, X4 has one more solid line than dotted line emanating 
from it, the conformal weight is unity. If we multiply back by the X24^x1^ factor removed 
previously, then we obtain an integral which is conformally invariant with respect to the 
external as well as internal points. 

The assumption of conformal invariance for the integrals immediately implies the "no- 
triangle" rule for the momentum-space diagrams. A loop with only three propagators would 
necessarily result in a negative weight for the x point corresponding to the loop momentum, 
because only three lines enter the dual diagram vertex. That negative weight can only be 
eliminated by additional denominator powers of x — that is, by additional propagators 
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FIG. 8: The rung-rule dual diagrams. A factor of st has been removed. 

which would turn the triangle subgraph into at least a box subgraph. 

By placing the dual diagrams on top of the original momentum-space diagrams, we can 
read off directly the change of variables between the Xij and the momenta: Xij is just the 
sum of the momenta of the each momentum-space line crossed by the dual line running from 
Xi to Xj. In the rung- insert ion rule, when a rung is inserted between two parallel lines with 
momenta h and I2, to go from an L-loop contribution to an {L + l)-loop one, the loops on 
either side of the parallel lines have acquired a new propagator. Hence each of their dual x 
vertices has a new solid line emanating from it. The rung-rule momentum-insertion factor 
of + ^2)^ is represented by a dotted line stretching between the two vertices, so it restores 
the conformal invariance for those two loops. This property may help to explain the form 
of the rung rule. 

Let us now focus on four-loop four-point integrals. We may find all conformally-invariant 
integrals by drawing the set of all dual diagrams that have conformal weight zero with respect 
to all Xi- Again to prevent cluttering the diagrams in figs. [8] and [9] with dotted lines we have 
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(dj) (f2) 

FIG. 9: The non-rung rule dual diagrams. A factor of st has been removed. 




(d') (d,) 

FIG. 10: The two four- loop dual integrals, in addition to those given in figs. [2] and [U that survive 
the requirement of conformal invariance. Both integrals are ruled out by two-particle cuts. In this 
case, no factor of st has been removed. 

removed a factor of st from the figures. The complete list of conformal four-loop integrals, as 
it turns out, contains the rung-rule diagrams of fig. [U the non-rung-rule diagrams of fig. [9l 
and the two extra integrals shown in fig. [TOl ^ 

Because conformal invariance in the sense discussed above implies the "no-triangle" rule, 
the complete list of candidate graphs for conformal integrals are given by the same no- 
triangle list described above, namely the diagrams in figs. [2], HI [S] and O Upon attempting 
to draw the dual diagrams to figs. [5]and[6l we find that in all cases except (ds), either it is 
impossible to add dotted or even solid lines so as to obtain a conformally-invariant integral, 
or the obtained conformally-invariant cases are equivalent to previously-included cases, or 

^ It is possible to dress diagram (f) in fig. [8] with dotted lines in a second way, but that dressing is related 
simply to the one shown by the exchange of legs ki and k^, or in other words s <-^- t. 
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FIG. 11: Two examples of dual diagrams which do not lead to new conformal integrals in the 
on-shell limit, for reasons discussed in the text. 

else conformal invariance can only be achieved by adding solid or dotted lines connecting 
neighboring external vertices. The latter lines are not admissible, however, because differ- 
ences of neighboring external x points correspond to individual external momenta fcj. The 
corresponding factors are then kf. Their presence would lead to an unwanted power-law 
vanishing or divergence of the integral in the on-shell limit. 

Figure [TT] illustrates two examples of dual diagrams that cannot be made conformal, or 
that reduce to previous cases. First consider the graph labeled (b4). The internal dual 
points all have weight four, so no dotted lines can attach to them. The dotted line shown 
between external points Xi and X3 reduces the Xi weight to zero and the X3 weight to one. It 
is nonvanishing in the massless limit. However, there is no other numerator factor available 
that is nonvanishing in this limit, to further reduce the conformal weights of external legs 
X3 and X4. In particular, = ^ in the on-shell limit. 

As a second example, consider the graph labeled (di) in fig. [Til Here there is one pentagon 
subgraph, and hence one internal x point to which a dotted line can attach. The dotted line 
shown can be used to reduce the conformal weight of X4 from three to two, which would then 
balance its weight with that of the opposite external point X2- (Balanced opposing weights 
can always be reduced to zero using powers of s = X24 or t = xl^.) However, this choice 
of dotted line merely cancels the propagator that it crosses, and thereby reduces the (di) 
graph to the (d2) graph in fig. [9], which we already know is conformally invariant and present 
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in the four-loop planar amplitude. Without using the dotted line shown, it is impossible 
to balance the X4 conformal weight, in the massless limit, so the only surviving possibility 
reduces to an existing integral. 

It is curious that the two conformally-invariant integrals represented in fig. [T0|, which are 
not present in the planar four-loop four-point amplitude, can be distinguished from the eight 
in figs. [S] and [H] that are present, by the fact that they do not have explicit overall factors 
of both s and t. As drawn in fig. [TUl they have three powers of s, but no powers of t. The 
integral (d') has the same basic topology as (d) in fig. [9l but the dotted lines emanating from 
the two pentagon loops are connected to external legs Xi and X3 instead of to each other. 
So it has two "rung-rule-type" numerator factors involving the squares of the sums of three 
loop momenta, but no power of t. At the moment, however, we have no good argument why 
explicit factors of both s and t have to be present in order that an integral be present in the 
four-point amplitude. 



VI. ANALYTIC AND NUMERICAL RESULTS 



Our next task is to evaluate the integrals entering the four-loop planar amplitude (14. ip 
in a Laurent expansion around e = 0. For the case at hand, massless gluon-gluon scattering 
in A/" = 4 super- Yang-Mills theory, all of the integrals encountered can be evaluated through 
three loops in terms of a class of functions known as harmonic polylogarithms (HPLs) [63 1. 
We expect this class of functions to continue to suffice at four loops. We know it suffices 
through 0{e~^), for which we have analytic results. 



A. Analytic expressions through 0{e 



The a naly tic results for the four-loop integrals were obtained with the help of the MB 
program 62]. We let x = —t/s and L = ln(— x). Through 0{e~'^) the results, expressed in 
terms of the HPLs defined in appendix [Bl are. 
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Using eq. (14. ip . together with the above results for the integrals, the total four- loop planar 
amplitude, m['^\ has the expansion. 
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In tableiniwe present the numerical values of the eight integrals appearing in the four-loop 
planar amplitude, through (9(e~^). The values through 0{t~'^) can be found easily from the 
analytic expression given above. The numerical values for e~'^ and were obtained using 
the CUBA numerical integration package 85|], which is incorporated into the MB program 62|. 
In the table, we also give the total value of the amplitude Mf\ according to eq. (14. ip . 
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TABLE II: Numerical values of individual four-loop integrals, and M| , at (s,t) = (—1, — !)• The 
uncertainties at orders and are indicated in parentheses. (The presence of two digits in 
parentheses signifies the uncertainty in the last two digits of the central value.) 



Integral 


e-8 


e-7 


e-6 


e-5 


e-4 


e-3 


e-2 


(a) 


4/9 

1 





-4.27225931 


-11.30789527 


-18.44325855 


-58.84504(10) 


-180.852(3) 


(b) 


4/9 





-4.82057067 


-10.53107909 


3.00827162 


67.62584(13) 


190.235(3) 


(c) 


4/9 





-5.30034310 


-10.38082198 


12.55376125 


81.91311(64) 


99.292(5) 


(d) 


4/9 





-3.86102580 


-7.70172456 


-16.94003184 


-80.03212(51) 


-52.555(21) 


(e) 


4/9 





-4.47787607 


-8.45252237 


-4.13769237 


-11.60392(20) 


28.823(9) 


(f) 


8/9 





-10.73776405 


-16.90406921 


46.68731190 


219.08111(12) 


364.167(7) 


(d2) 











-0.80137127 


2.48032408 


36.23672(11) 


132.811(13) 


(f2) 


16/9 





-21.47552810 


-35.41088096 


92.92365579 


521.48787(31) 


1314.856(12) 




2/3 





-7.12804762 


-13.64293336 


2.64276920 


27.34123(13) 


33.278(7) 



B. Values of lower-loop amplitudes at (s,t) = (—1, —1) 

Next we need to compare our results for eq. (14. ip with the prediction (12.251) based on the 
known structure of the infrared poles. To do this, we evaluate the lower-loop amplitudes, 
using formulas from ref. Q], at the symmetric kinematical point (s,t) = (—1, —1), through 
the accuracy needed to evaluate eq. (12.251) to C(e^^). At this point, a limited number 
of analytic expressions appear, built out of In 2, tt, Cs; Li4(|), Li5(|), Li6(|), and the 



harmonic sum 



86|, 



S6 



°° ( — l)*i *^ ( — 1)'2 

5({-5, -1}, oo) = ^ XI = 0.98744142640329971377 .... (6.10) 



il=l ^ 12 = 1 

r(i)/ 



The one-loop amplitude Mi '{s,t]e) in eq. (13. ip . evaluated at {s,t) = (—1,-1), has the 
e-expansion, 
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-4 + 6.5797362673929057461 

+ 11.133742693869288271 e + 7.1556624851455749140 

- 5.760188577405266135 - 23.794568007684383734 + C(e^) . 

(6.12) 

r(2) 



The two-loop amplitude {s,t; e) in eq. 03.31) . is given by 
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V ^^2^ 30 9 360 12 80^^ 

+ e2 ^ 4 Li ( 1 ) + Li (1) + ^ 1^6 2 + ![! 2 - ^tt^ In^ 2 

V 6V2^^ 3 ^^2^^180 ^ 4 240 

307 2. , 4319,2 541 .\ 

+ vrVg n 2 + Cl + O 6.13 

24 288 6480 / ^ ^ ^ ' 

2 12.337005501361698274 



t 

+ 30.635795346547106621 e + 68.218654436238118625 + ^(e^) . 

(6.14) 

The three-loop amplitude M\ {s,t; e) in eq. (13. 5p is given by 



Mf (-l,-l:e) = -— + 



4 77r 



+ _L(_4Li4)-iln^2 + ^vr2ln^2 



89 4 

TT 



1/ . /is 1 , 2 2, s 73 4, 3779 2. 8621 \ 

+ - 4 Lis i In^ 2 + -vr^ In^ 2 vr^ In 2 TrVg + C5 

eV ^^2^ 3Q g 3gQ 432 120^/ 



+ 14s6 - 4 Li6(i) - -7r2 Li4(i) - — In^ 2 - — 7r2 In^ 2 + — tt^ In^ 2 
^ ^^2^ g 4V2; ^gQ 3gQ 



-23 7r2C3ln2 - + ^^n' + 0{e) (6.15) 

S3 ^44 S3 -r 233280 ^ ^ ^ ^ ^ 

4 11.514538467937585056 
+ 



3e6 

21.065428484578982255 1.6228781926783846589 
+ 7^ 7 
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40.219043687209842734 , 

67.305777557207060997+ C(e) . (6.16) 

e 

The iterative formula for the four-loop amplitude in terms of the lower-loop amplitudes is 
given in eq. (12.251) . Inserting the values at (s,t) = (—1,-1) of M^^ and we 

obtain 

2 13 TT^ 1 / 2 „ 68 



Mf (-l,-l;e) 



— - + ^( — TT^ ln2 - —6 

3 18 e6 e5 V 3 9^^ 

g4 Vs 4^2^ ^ 9 9 2592 

+ \ (J_ Li5(i) + 1 In^ 2 - -TT^ ln=^ 2 + — vr^ In 2 
e3 V 3 ^^2^ 45 27 135 



89 4 

-TT 



3 45 27 135 

251 2. 7469 

TT Cs 

36 120 



+ 1 f _i4sg + ! Li6(i) + — vr^ Li4(i) + — In^ 2 + — vr^ In^ 2 
^ 3 ^^2^ 9 ^^2^ 216 
607 4, 2 . 791 ^ 895^2 20759 g I ^u) 
-y \v? 2 + — vr^Ca In 2 + — C| - 77— -r^' - 0/0 



1080 36 432''' 102060 

+ O(e-i) (6.17) 
0.66666666666666666667 7.1280476230089812249 



g8 g6 



13.642933355332790075 2.6427691992903098962 



g5 ^4 



27.341074205440151100 

+ ^ 



+ ^(^29.611139840724282137- \fo^) + C'le"^) • (6.18) 

Comparing eq. ( 16. 18^ with the last row of table [Tll we verify the pole behavior of the 
four-loop amplitude m|'*''(— 1, — 1; e) precisely through 0{e~^). The agreement at 0{e~^) is 
good to 5 digits. At (9(e~^), we can extract the value of fl^\ We obtain, 

= -29.335 ± 0.052 . (6.19) 

This value should be compared with that predicted by Eden and Staudacher, from eq. (11.61) . 



JO 



73 

-n^ + = -26.404825523390660965 .... (6.20) 



Es 2520 

The results do not agree. The difference can be expressed as, 



/o^'^ = /o^'^ + Mr , (6.21) 

ES 

A/(S^^ = -2.930 ± 0.052 .... (6.22) 



.(4) 
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We can also parametrize the difference A/g as a multiple of the weight-6 expression 
Making this parametrization, we find that 

A/(^)=rC|, (6.23) 
r = -2.028 ±0.036... . (6.24) 

This result is quite suggestive. To about 1.5% precision on the value of the correction term 
A/q'*'', it is equal to a result which would have the net effect of flipping the sign of 

the Cl term in the Eden-Staudacher prediction while leaving the vr^ term unaltered. 

Of course, the ES prediction follows directly from an integral equation, and so flipping the 
sign of the Q term is not possible without other modifications. In section IVIIII we discuss 
possible reasons for the discrepancy. 



C. Cross checks at asymmetrical kinematical points 



In order to cross check our numerical evaluation of the integrals at the symmetric kine- 
matical point (s,t) = (—1,-1), as well as check the behavior of the (9(e~^) and (9(e~^) 
terms in the planar four-loop amplitude as a function of the scattering angle, we have per- 
formed the numerical analysis of the last subsection at three additional kinematical points, 
(s, t) = (-1, -2), (-1, -3) and (-1, -15). 

We have used the expressions for the lower-loop amplitudes in ref. [S] to numerically 
evaluate the infrared-based iterative formula fl2.25p at the asymmetric kinematical points 
(s, t) = (-1,-2), (-1,-3) and (-1,-15). The results are, 

0.92419624075 6.64759460909 4.23222757233 



Mf)(-l,-2;e) 



iter. 



2 

3^ 



15.89245103368 16.11914613046 
+ i + § 

+ ^(^1.31283053842- ^/^^^) + 0{e-^) 



(6.25) 



Mf (-l,-3;6) 



2 1.46481638489 

3^8 



5.92109866220 0.72092946726 
s + 



19.05722201166 4.86152575608 
+ i + § 

+ ^(-5.61581265989- ^/o^'^) + 0{e-^) 



(6.26) 
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Mf (-1,-15; 



2 3.61073360147 0.20548826868 14.13192416428 



+ s + 



iter. 3e8 e7 e6 

7.41700629511 48.55010675803 
+ -i ^3 

+ -^(43.61197714-^/0^')) + O(e-i) . (6.27) 

Numerical evaluation of the eight four-loop integrals from figs. [2] and H] gives the fol- 
lowing results for the amplitude m[^^ (we omit the 1/e® through 1/e^ poles, as they agree 
analytically) , 



Mr'(-l,-2;£) = 0(e'»...£-') + — — + / +0{e-'), (6,28) 

Mr)(-l,-3;.) = 0(.-«. . ..-^) + ^-8617 ±0.0003 _ 1.943 ±oW ^ ^ ^ 

f4W ^ 8 4^ 48.5499 ± 0.0002 47.29 ±0.02 
MfVl,-15;e = O e-'-.-e"' + + • 6.30 

Comparing these sets of numbers at 0{e^^), we observe good agreement at all points; the 
{s,t) = (—1, —2) point is slightly off, at 1.9 cr, but the other two points are within la. 

At (9(e~^), we can express the agreement in terms of the parameter r introduced in 
eq. (16.231) . At the asymmetric kinematical points, we extract the values, 

r = -2.059 ±0.036, (s, t) = (-1, -2), (6.31) 
r = -2.062 ±0.045, (s, t) = (-1, -3), (6.32) 
r = -2.074 ±0.104, (s, t) = (-1, -15). (6.33) 

These values are all consistent, within errors, with the value fl6.24p extracted at {s,t) = 
(-1,-1). (The values at different points, however, have an unknown correlation between 
them, because the integrals contain pieces that are independent of the kinematics, and the 
numerical integration for each value of (s, t) was performed with the same sequence of quasi- 
random integration points. This means the results for r from the various kinematic points 
cannot be combined to reduce the error.) 

In summary, our numerical integration of the four-loop planar integrand results in a value 
of the cusp anomalous dimension, 

/o(a) = a-—a^ + —7T^a^-(^n''-{l + r)Ci)a^+ (6.34) 
■'^^ ' 6 180 V2520 ^ '^^) ^ ' 

where r is given in eqs. (16.241) . (I6.3ip . (16.321) and (I6.33p . All values are consistent with the 
appealing value of r = —2, which corresponds to the value /? = ^3 for the dressing-factor 
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parameter /3 in eq. fll.7p . As noted above, this value would merely flip the sign of the 
el term in the ES prediction fll.Sp . However, we obviously cannot exclude, on numerical 
grounds alone, nearby rational or transcendental numbers. For example, it is conceivable 
that r takes on the value r = —5/(2(^3) = —2.0797. ... This value would correspond to 
a rational dressing- factor parameter, (3 = 5/4, and would violate the KLOV maximum- 
transcendent ality principle. On the other hand, additional evidence points toward r = —2 
as the correct analytical value, as we shall discuss in the next section. 



VII. ESTIMATING STRONG-COUPLING BEHAVIOR 



Kotikov, Lipatov and Velizhanin (KLV) 34] made an intriguing proposal for approxi 



mating the cusp anomalous dimension (or equivalently, fo{a)), for all values of the coupling 
a. They suggested combining perturbative information with the knowledge from string the- 
ory [69] that at large values of the coupling /q has square- root behavior, /q ~ Vd,. They 
proposed the following approximate relation as a means for incorporating the known analytic 
behavior, 

2n 

a" = ^a[/o(a)r, (7.1) 

r=n 

where the constants Cr can be fixed using perturbative information. As we shall discuss 
below, they can also be fixed using strong-coupling information. As more information be- 
comes available the integer n can be increased. The strong-coupling square-root behavior 
of /o is automatically imposed by the fact that a" ~ C2n[/o(d)]^"' at large a. Similarly, the 
weak-coupling linear behavior follows from d" ~ Cn]/o(Q')]" small a. 

KLV used the approximation (17. ip for n = 1, together with the one- and two-loop ex- 
pressions for the cusp anomalous dimension, to write (for the case of a supersymmetric 
regulator) 

a = /o + ^ (for ■ (7.2) 

6 

This formula makes the weak-coupling predictions (beyond two loops), 

f, = a--a' + -d'- — vr^ a" + — vr^ d' - -^vr^o a' + . . . , 
•' 6 18 216 648 1296 

as d ^ 0, (7.3) 
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and it predicts coefficients in the strong-coupling expansion, as 

2^3 /a 3 



fo = \l^ - ^ + 0{d-'/') (7.4) 



TT V 2 TT 

1.1027 - 0.30396 + ^(a^^/^) . (7.5) 



The coefficients of the leading [33|, l69|, l70|] and subleading [7l|] terms in this expansion are 



predicted from string theory to be, 

^ - 0.16547670011448 + ^(a-^/^) . (7.7) 

As noted by KLV, the leading coefficient is estimated correctly to 10% by the formula (17.21) . 
The subleading coefficient is off by almost a factor of two, however. 

What happens as we incorporate more perturbative information? Using the three-loop 
value for the cusp anomalous dimension, and setting n = 2 in eq. (17.11) ). gives the approxi- 
mation 

= (/o)^ + y (/of + ^ (/or ■ (7.8) 
This approximation makes the weak-coupling predictions (beyond three loops), 

f ^ a2 I JJl^4 -3 _6 -4 I 329 g 5 169 iQ -6 , 

to = a a H Tc a tt a H vr a vr a + . . . , 

6 180 1080 21600 19440 

as a ^ 0, (7.9) 
and has the strong-coupling expansion, 

/o = -15^/^\/f - + 0{a''/') (7.10) 

TT V 2 TT^ 

1.2529 - 0.50661 + 0{a-^/'^) . (7.11) 

For both the leading and next-to-leading coefficients in the strong-coupling expansion, the 
three-loop approximation (17.81) leads to a larger disagreement with the string prediction (17. 7p 
than does the two- loop version (17. 2p . We also note that the numerical value of the four-loop 
coefficient predicted by eq. (17. 8p is —27.595, which is a bit closer to our result than is the 
ES prediction, but still about 6% off. 
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Despite the somewhat discouraging results from including the three-loop values, we pro- 
ceed to incorporate our four-loop cusp anomalous dimension into the n = 3 version of the 
approximation, obtaining 

a' = ifo)' + y ifo)' + ^ ifof + - 3(1 + r)C|) (/o)^ • (7.12) 

Here we have introduced the same coefficient r defined in eq. (16.231) . which is constrained 
to be quite close to —2 by our numerical result fl6.24p . The weak-coupling expansion of this 
formula predicts (beyond four loops), 

/o = a - - a2 + —7T^ a^-( ^vr^ - (1 + r)Ci] 
6 180 V2520 ^ 

+ TT^ - -(1 + r)7rV.? a - (1 + r)TT^d ] + . . . , 

V113400 3^ ' V453600 10^ ^ J J ^ 

as a 0. (7.13) 
The strong-coupling expansion is given by, 

/o = - ^ + — (^ + {1 + rxi) a-''/' ./f + 0{a-') , (7.14) 

•'^ V 2 720a 256 V2835 ^ \ a ^ ^ ' ^ ' 

where 

a = --(l + r)d, (7.15) 

3024 8^ ^ ^ 

and we have given one more term in the expansion than before. Curiously, with the ES 
value r = 0, the approximate relation fl7.12p breaks down, because a is negative, and hence 
Q,-i/6 jg j^Q^ real. 

For r = — 2, formula (17.141) becomes, 

/o ^ 1.02550 yi - 0.157356 - 0.0562398 + 0{d-^) . (7.16) 

The numerical agreement between eq. (17.161) and the string-theory result in eq. (17. 7p is quite 
impressive: The leading coefficient agrees within 2.6%, and the subleading coefficient within 
5%. The coefficient of the term proportional to \/2/d in eq. (17.16^ is fairly small. As we 
discuss below, an improved estimate suggests that it may be considerably smaller, or perhaps 
even vanish. 

In the KLV type of approximation, the predicted value of the strong-coupling coeffi- 
cients depends quite sensitively on the value of the four-loop contribution to the anomalous 
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dimension. For example, if we scale the numerical value of the four-loop contribution as 
follows, 

instead of eq. fl7.16p . we find at strong coupling, 

/„ . (0.8597725 + 10.9855 -^^5^3^^^^ + O(a-V^) , (7.18) 

which exhibits a strong sensitivity under just a few percent change in the four-loop contri- 
bution. More generally, the sensitivity of the strong-coupling prediction to the higher-loop 
orders used in a KLV approximation, allows us to test whether a given ansatz appears 
compatible with strong coupling, as we discuss below. 

In fig. [12] we plot these estimates as a function of the coupling, and we also display 
the strong-coupling limit (17.71) predicted by string theory. As noted above, the approxima- 
tion (17. 2p using only two-loop information works quite well, in fact better than the three-loop 
approximation (17.81) . However, the behavior of the four-loop formula (17.121) is clearly ex- 
tremely close to the string theory prediction. For this curve, or rather band, we have varied 
the parameter r between —2.05 and —1.95, consistent with eq. (I6.24p . 

We have constructed two further approximations by incorporating the knowledge of the 
precise strong-coupling coefficients. Matching the four-loop perturbative information and 
the leading strong-coupling coefficient gives the approximation 



a' = {for + + ar + (^^' - 4(1 + oci ) (hy + le (/, 



0) ■ 

(7.19) 



It has the weak-coupling expansion, 

6 180 V2520 ^ 

.2,2 , \ ^5 



/ 4747 s 3,^ , A 

V302400 2^ ^ J 

f 5023 10 19,^ , 20 2\ .6 

- 7r^° (1 + r)7rV| vr^ \ + . . . 

V 544320 12^ ^ 3 y 



as a -> 0. (7.20) 

If we add the next-to- leading strong-coupling coefficient as another constraint, we obtain 
the approximation, 

l5 



a 



Uof + {for + y {for + (^^^ - 5(1 + r)C|) {foY 
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Approximate Cusp Anomalous Dimension in Planar MSYM 
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FIG. 12: Approximations to the cusp anomalous dimension in planar MSYM based on the for- 
mula (j7.ip and perturbative information through two loops (dotted line), three loops (dot-dash 
line) and four loops (solid line). Also shown is the strong-coupling prediction from string theory 
(dashed line). 
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with the weak-coupling expansion, 



Jo = a a + -77^ TT a 
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vr« - (1 + r)Q a' 
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TT 



as a — i> 0. 



(7.21) 



(7.22) 



We may also use the improved approximation (17.2 ip to determine subleading terms in 
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the strong-coupling expansion. We have 

~ /a 31n2 / 17 81 , 2 1 . ^ '?\ u . X 

fo = \ ( 7^ ^ il + r)d]\ - + 0(a-^). 7.23) 

•^^ V 2 477 V161280 327r2 64^ '^^J\d ^ ^ ^ ^ ^ 

For r = — 2, this expression evaluates to 

^ . 12 



/o ~ \/^ - 0.16548 - 0.000693 + C(a"^) . (7.24) 



The first two numerical coefficients automatically reproduce the input [33|, |69|, |70|, lDJ from 
string theory, eq. (17. 7p . so the content of this equation is in the coefficient of ^/2/a. Quite 
interestingly, this coefficient is very small, suggesting that it may even vanish in the exact 
expression. It is an intriguing question whether such a tiny or vanishing result could be 
obtained from string theory. 

We plot these two approximations, for r = —2, along with the four-loop approxima- 
tion (17.121) . in fig. [131 Obviously, they are extremely close to one another. As a measure of 
that, the values they predict for the five-loop planar cusp anomalous dimension (for r = —2) 
are 167.03, 166.34 and 165.25, corresponding to the three approximations in eqs. (17.131) . 
(I7.20p and (I7.22p . We expect the last of these numbers to be the most accurate, given that 
eq. (17.22P uses all available information at both strong and weak coupling. 

Can we use these approximations to guide corrections to the ES prediction at higher 
loops? As a ffist step, consider the five-loop ES prediction in eq. (II. 6p . The ES numerical 
value is 131.22, which disagrees significantly with the prediction of our approximate formula 
(]7.22p . However, we can modify the five- loop ES prediction (11.60 in a manner analogous to 
the modification required at four loops to fit the value r = —2. We flip the signs of the 
terms containing C„ values for odd n, and leave untouched the terms containing only even ( 
values (those terms containing only 7r's and rational numbers). Then the terms containing 
(^3 and/or in the five-loop coefficient in eq. (II. 6p acquire the same sign as the 7r^ term, 
instead of the opposite sign, giving an analytic expression through five loops: 

^modificd/-\ - -2 , -'--'-4-3 f '^'^ ^6 , 

In (a) = a a H n a — n + C2 ]a 

■'0 ^ ' Q 180 V252O 

+ f J!L7r8 + -C| + 10C3C5V5+ (7.25) 
V 56700 3 J ^ ' 

The numerical value of the five-loop coefficient in eq. (I7.25p . 165.65, agrees with the five- 
loop coefficient in our approximate formula ( 17.22P to a remarkable 0.25%. This excellent 
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Approximate Cusp Anomalous Dimension in Planar MSYM 
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FIG. 13: Approximations to the cusp anomalous dimension in planar MSYM based on the for- 
mula (j7.ip . four-loop perturbative information, and zero, one or two constants from the strong- 
coupling expansion. 



agreement in turn reinforces the notion that r = —2 gives the correct analytical value 
of the four- loop coefficient. (We remark that a similar procedure at three loops, using the 
approximation (17. ip for n = 4 and incorporating the two leading strong-coupling coefficients, 
works very well to estimate the next perturbative term: It predicts a four-loop coefficient of 
30.22, which compares nicely with our result (I6.19p .) 

We now continue this procedure to higher loops, comparing various predictions to our 
approximate formula. To obtain the ES prediction to higher-loop order, we use the integral 
equation from which eq. (II. 6p is derived |31| . 

Ji(V2at) 



/(a) = a- Aa^ I dta{t) 
'o 



'2at 



(7.26) 
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where the fluctuation density ait) is obtained by solving the integral equation, 



t 

a{t) 



e* - 1 
with the kernel, 



Ji(v^t) 



POO 

Jo 



(7.27) 



K^,^,^ = JJM^^mm^ (7.28) 

where Jo and Ji are standard Bessel functions. 

Solving this equation perturbatively through 12 loops, we find the numerical values of 
the ES prediction to be, 

/ES = a - 1.6449 + 5.9528 a^- 26.405 + 131.22 a^- 705.54 a*^ 
+ 4021.9 - 23974. a^ + 1.4800 10^ a^ - 9.3958 10^ a^° 
+ 6.1024 10^ a" -4.0387 10^ a^^ + -- - , 

as a ^ 0. (7.29) 

The values in eq. fl7.29p may be contrasted to the ones obtained from the weak coupling 
expansion of our approximate formula (17.211) (with r = —2), 

y.approx ^ . _ ^ g^^g .2 ^ 5 9528 - 29.295 a^ + 165.25 - 1002.7 a^ 
+ 6379.3 a^ - 41997. a® + 2.8371 10^ a^ - 1.9555 10^ a^° 
+ 1.3699 10^ a" - 9.7237 10^ a^^ + ■ ■ ■ , 

as a ^ 0. (7.30) 

The agreement between eqs. (I7.29P and (17.301) at five loops and beyond is rather poor. This 
disagreement is not surprising, given that the ES four-loop value, used as input into our 
approximate formula (I7.30p . differs significantly from our calculation. 

We have evaluated further terms in the weak-coupling expansion of our approximate 
formula (I7.2ip . up to 75 loops. The ratio of the n^^ term to the {n — term in the series 
slowly settles down to a value near (—8) — at 75 loops it is (—7.95) — suggesting a radius 
of convergence of 1/8, and a nearest singularity on the negative real axis at dc = —1/8. This 
value does appear to agree with the location of the nearest singularity in the original ES 



equation [471, l87 | 



What happens if we generalize the four- and five-loop sign flips in the ES prediction, and 
require that all contributions at a given loop order come in with the same sign? Doing so, 
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we obtain, 

ymodified ES ^ ^ _ I Q4^^Q ^2 ^ 5 9528 - 29.295 CL^ + 165.65 - 1007.2 
+ 6404.7 - 42020. a® + 2.8223 10^ - 1.9307 10*^ 
+ 1.3406 10^ a" - 9.4226 10^ a^^ + ■ ■ • , 

as a — 0. (7-31) 

This modified formula is much closer to our approximate expression (17.301) . differing even at 
12 loops by only 3 percent. Of course, we have no reason to expect this naive modification 
of signs in the ES formula to be correct to all loop orders, though it appears to get the 
numerically largest contributions correct. Indeed, if we use this series to systematically 
construct KLV approximations with larger values of n in eq. (17.11) as more terms are kept, 
the large a coefficients do not settle to the string values (17.61) . For example, using the weak- 
coupling series (I7.3ip up to eight loops as input to the KLV approximation (17.11) with n = 7 
leads to a prediction at strong coupling, 

a . .^^^-1/2^ 



fo ^ 0.95968 -\/- - 0.081182 + 0{d-''^) . (7.32) 



t is noteworthy that the series (I7.3ip corresponds to the contemporaneous proposal 
. 47I. (We have checked agreement through 30 loops.) In order to improve the 



which compares poorly with either the string result (17.71) or with the four-loop prediction 

dnsi). 

in ref. 

agreement of the higher-loop KLV approximations with the string theory strong-coupling 
coefficients, further modifications to the proposal are needed. Similar conclusions follow 
from the investigation of a sequence of Fade approximants, as discussed below. 

The surprisingly good agreement of eqs. (17.301) and (17.311) does, however, suggest that a 
simple repair of the integral equation is possible. As a trivial example, by modifying the 
kernel in eq. (I7.27P to K{\/2d,t, —\/2d,t'), we obtain 

ymodified K ^ ^ _ ^ g^^g ^2 ^ 5 9528 - 29.295 + 165.65 - 1011.9 
+ 6490.0 - 43050. a® + 2.9271 10^ - 2.0282 10^ 
+ 1.4265 10^ a" - 1.0156 10® a^^ + ■ ■ ■ , 

as a — 0. (7.33) 

which is in fairly good agreement with our approximate result. Although this ad hoc modifi- 
cation should only be taken as an illustration, it does show how one can use our approximate 
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formula fl7.30p to guide corrections to the ES integral equation. (It also suffers from the 
problem that strong-coupling extrapolations do not match the string values.) A proposed 
modification of the integral equation, valid through (9(a^) and leading to a modification in 



the anomalous dimension at 0{a^), is given in eqs. (89) and (91) of ref. 3l|]. The choice 
f3 = corresponds to r = —2 in eq. (16.23^ . as can be seen from eq. (92) of the same refer- 
ence. This is also the choice preferred by crossing symmetry in the strong-coupling limit of 



the dressing factor, within the family of modifications newly studied in ref. j47l ]. 

It is useful to investigate other approximation schemes, to make sure that our results 
are not significantly biased by the form of the KLV approximate formula (17.11) . A well- 
known method for incorporating information from different expansion regions — here from 
both weak and strong coupling — is that of Pade approximants, which fit a function to 
a ratio of polynomials in a suitable variable. In the complex a plane, we expect that 
fo{a) has poles or branch cuts. As discussed earlier, there is evidence from the behavior 
of the KLV approximation and the ES kernel that the singularity nearest the origin is at 
cic = —1/8. For the purposes of constructing the Pade approximant, however, we will model 
the singularities by a branch cut terminating on the negative real axis at the point —^"^/S. 
This assumption leads us to introduce a variable u = a/1 + 8a/^^, and define the [(m+l)/m] 
Pade approximant by 

ft'-''^-\u) = {u- ,^No + mn + ... + N.^u- 
° V ; V ^ 1 + Diu + ... + DmU"" ^ ' 

The relative degrees of the numerator and denominator polynomials are fixed by the strong- 
coupling requirement /o ~ \fh ~ u. The factor of {u — 1) comes from the vanishing of /o at 
a = 0. The remaining 2m + 1 constants can be fixed by perturbative data, or by a mix of 
weak- and strong-coupling data. Note that when m increases by 1, two more input numbers 
are required. 

Consider first the [3/2] Pade approximant, 

(7.35) 

which is implicitly a function of ^, which parametrizes the location of the cut termination 
on the negative real axis. We solve for A^o? A^i, A'"2, D\ and D2 as a function of ^, using the 
two strong-coupling coefficients in eq. (17. 6p . plus the one-, two- and three-loop coefficients. 
We then expand the result in a in order to estimate the four-loop coefficient. This procedure 
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[3/2] Fade Prediction of Four— Loop Coefficient 

—26 I 1 1 1 1 1 1 1 1 1 1 — I I I I 1 1 1 1 1 1 1 1 1 1 1 r 




I I I I I I I I Li I I I I I I / I I I I I I I I I I I I 

0.0 0.5 1.0 1.5 2.0 2.5 

FIG. 14: Value estimated for the four-loop cusp anomalous dimension /q^^ by the [3/2] Fade 
approximant (j7.35p . as a function of the parameter ^ controlling the point at which the cut is 
assumed to terminate, a = — ^^/8. 

will give us information about ^. The four- loop coefficient is plotted in fig. O The value 
^ = 1 is picked out rather clearly by this plot, as the (approximate) first location where 
the curve crosses the value ff^ = -(73 7rV2520 + C|) = -29.2947 (for r = -2). At ^ = 1, 
formula fl7.35p predicts —29.521. The second crossing is a spurious feature of the Fade, and 
the third crossing, near ^ = 2, is highly implausible, based on the strong growth in the 
perturbative coefficients already at four loops, not to mention at higher orders in the KLV 
approximate formulas. 

In fig. [15] we plot the ratio of the [3/2] Fade to the KLV approximate formula (17.211) . 
for various values of ^ near 1 that produce four-loop coefficients that are not too far from 
—29.2947. Except for the pathological case of ^ = 1.05, all the ratios are within about 1% 
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of unity. For the preferred value of = 1, the [3/2] Fade agrees with the KLV approximate 
formula everywhere to within about 0.1%. We also investigated the [4/3] Fade, for ^ = 1. 
This approximant requires two more input terms, namely the four- and five-loop coefficients. 
We used the "sign-flipped" values in eq. (17.251) . The resulting expression predicts the six-loop 
coefficient to be —1005.5, in quite good agreement with the value of —1002.7 in eq. (I7.30p . On 
other hand, a plot of this function reveals a singularity on the positive real axis, at u = 4.791, 
which is an artifact of a positive root to the cubic denominator polynomial. (The quadratic 
denominator polynomial in the [3/2] Fade for ^ = 1 has complex roots, relatively far from 
the real axis.) We also constructed a sequence of Fade approximants (I7.34p to see how the 
strong-coupling coefficients are predicted by the proposed sign-flipped sequence (I7.25p . As 
with the KLV approximation, we do not flnd any convergence to the string values ( 17. 6p . For 
example, the [4/3] Fade, based on the flrst seven loops in eq. (I7.25p . estimates 

/o ^ 0.95696 - 0.061164 + 0{d^^/^) , (7.36) 

which is close to the KLV estimation (17.321) . and signiflcantly different from the string values. 

In summary, the different nature, yet very similar numerical results from the Fade approx- 
imation method gives us confidence that the KLV approximation (I7.2ip is good to better 
than a percent for all values of the coupling. It also gives confirming evidence that the 
singularity nearest the origin is located at = —1/8. 



VIII. ANALYSIS AND CONCLUSIONS 



Using unitarity, we constructed the planar four-loop four-gluon amplitude in MSYM, and 
verified the structure of its infrared divergences through 1/e^. At order 1/e^ we were able 
to numerically test the prediction of Eden and Staudacher for the four-loop cusp anomalous 
dimension. Our result disagrees with their predictions. Using an approximate interpolating 
formula due to Kotikov, Lipatov and Velizhanin 3J], with the first four loop orders as input, 
we are able to estimate the first two coefficients of the strong-coupling expansion, coming 
to within 2.6 and 5 percent of the string theory prediction. Improving the approximate 
formula by using the string predictions as additional input, we obtain a formula for the 
cusp anomalous dimension which we expect to be valid for all couplings, to better than one 
percent. We have confirmed this using Fade approximations. Not surprisingly, given the 
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FIG. 15: Ratio of the [3/2] Fade approximant (I7.35p . as a function of the parameter ^, to the KLV 
approximate formula (|7.21|) . as a function of the coupUng a. 



disagreement at four loops, at higher orders the weak-coupling expansion of our approximate 
formula fl7.2ip has large disagreements with the ES prediction (17.291) . 

There are a number of conceivable reasons for the discrepancy at four loops, which we 
shall examine briefly in turn: 

• Different definition of the coupling g. 

• Difference between the cusp anomalous dimension defined for space-like vs. time-like 
kinematics. 

• Breakdown of integrability. 



The wrapping problem. 
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• The dressing factor. 

We expect the latter reasons to be more hkely than the former. 

In principle our result could differ from the ES result simply because we are expanding 
in coupling constants that differ beginning at three loops, say a! = a + xqQ for the 
appropriate value of xq. We use the classical coupling as our expansion parameter, so if 
this explanation is correct, the asymptotic Bethe ansatz would have to be using a different 



coupling. For example, in the magnon dispersion relation, discussed in eq. (2.18) of ref. [67|, 
a parameter is adjusted to match the perturbative coupling of MSYM, via the relation 
a(3 = a/2. (Here a and P describe the action of two central charges adjoined to an SU(2|2) 
algebra acting on the spin chain.) If this simple substitution were to have corrections 
starting at order a^, it would lead to a four-loop discrepancy between the ES prediction 
from integrability and our direct calculation in MSYM. 

Suppose our result and the ES result corresponded to cusp anomalous dimensions of 
different signature, i.e. one space-like and one time-like. Could this difference account for 



the different values? According to the results of Dokhshitzer, Marchesini and Salam |88| . 
it cannot, because the leading large-x limits of the space-like and time-like DGLAP kernels 
should be identical, to all loop orders. 

Integrability of the dilatation operator, interpreted as a spin model Hamiltonian, has not 
been proven to hold to all orders. However, the fact that there are very similar integrable 
structures at very strong gauge coupling — in the classical sigma model {25! and even in its 



quantum corrections [26|] — suggests that integrability should persist. 

Generically, a wrapping problem can occur if one tries to apply a Bethe ansatz solution 
for a spin-interaction that has longer range than the number of spin sites, which is equal 
to the twist J. Because the range of the interaction increases with the loop order, at fixed 
J this problem becomes more severe with increasing loop order. It may well be that the 
wrapping problem prevents the asymptotic Bethe ansatz from being applied to short (twist- 
two) operators, even if it can be used for longer ones. The wrapping problem is generically 
supposed to arise at order a"^"^ for operators of twist J. It is only other symmetries that 
protect the one-, two- and three-loop cusp anomalous dimensions from being affected by the 
wrapping problem. 

It is quite possible that the discrepancy is resolved by the so-called "dressing factor", 
an overall phase for the AdS/CFT ^-matrix, which is consistent with integrability, the 
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PSU(2,2|4) symmetry, crossing symmetry, etc. The strong-coupling expansion of the dress- 
ing factor is known to be nontrivial. Indeed, the first two terms in the semi-classical expan- 
sion of the dressing factor on the string side have been worked out [4^ 



44| . There have 
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also been very interesting recent analyses of the properties of the dressing factor under the 
worldsheet crossing symmetry 45|, |46l]. However, it remains unclear at which order in the 
weak-gauge-coupling expansion the dressing factor becomes non-trivial. If the reason for the 
discrepancy is indeed the dressing factor, then we now know it begins at four loops. 

Eden and Staudacher made a specific proposal for modifying the asymptotic Bethe ansatz. 
They proposed a dressing factor containing a parameter j3, thereby modifying the kernel of 
their integral equation. This leads to a shift in the four-loop cusp anomalous dimension 
proportional to /3. If this proposal is correct, we determine the value of this parameter to 
be close to, if not exactly, (3 = (3. 

Anomalous dimensions of many other classes of operators are linked through the 
PSU(2,2|4) symmetry, and are therefore affected by a nontrivial dressing factor. As just 
one example, the anomalous dimension of the operator O = Tr(X^Z^) + . . ., where X and 
Z are two of the three complex scalar fields in MSYM, is altered by an amount proportional 
to P, to have the form |3l| . 

■115 



(8.1) 



Thus the Eden-Staudacher proposal for the resolution of the discrepancy could be tested 
by a direct computation of the four-loop anomalous dimension of O.^ For /? = ^3, it would 
appear to result in a non-uniform transcendentality for this quantity [31]. (A caveat is that 
one needs to account also for the transcendentality assignment of harmonic sums (35], which 
are sums of rational numbers.) As pointed out by Eden and Staudacher, a nonzero value of 
/?, through its effect on the anomalous dimension of O, apparently rules out the Hubbard 
model Hamiltonian as a candidate for the SU (2) dilatation operator beyond three loops 89|. 

The Eden-Staudacher modification is the lowest-order term of a general form for the 
dressing factor 5*12 = exp(2i6'(a:i, X2)), with 



00 00 



{Xk)qr{Xj)) 



(8.2) 



r=2 v=r-\-l 

suggested in the literature Q, Q, 165|, 



66, 



6:^, 



681]. In this equation, the qr are the 
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spin-chain charges, 

1 



and the xf are rapidities entering into the Bethe ansatz by which the spin-chain has been 
solved at lower orders. 

All known quantities, anomalous dimensions or amplitudes, in the Af = 4 supersymmetric 
theory, have a uniform transcendentality. That is, all polylogarithms or zeta constants at a 
given order in the perturbative expansion have the same polylog weight or transcendentality. 
The ES proposal deviates from this observed property in some quantities, such as eq. (18.11) . 

Is the dressing phase in eq. fl8.2l) general enough to accommodate our expectations for the 
A/" = 4 theory? As noted by Eden and Staudacher, maintaining uniform transcendentality 
with purely rational coefficients requires every power of t or t' in their integral equation to 
come along with a power of the coupling g. Now, the coefficient of the C| terms, which 
we wish to modify, arises from terms in the kernel proportional to tt'. This term is odd 
in t'. Such a term cannot arise from Fourier transforming a lone charge q, because we 
are interested in symmetric densities, and hence will ultimately take the symmetric part 
of this transform. The symmetric part is either zero (for with r odd) or a symmetric 
function of t'. In eq. (18. 2p . every term is linear in a charge q{xj), and so it cannot lead to 
the modifications we need. We must therefore seek a more general dressing factor, or else 
introduce non-rational coefficients. 



Very interestingly, in a contemporaneous paper, Beisert, Eden and Staudacher j47|] have 
done just the latter, in such a way that, at least for the anomalous dimension under consid- 
eration, uniform transcendentality is maintained. This leads to a modified integral equation 
of the ES type for the cusp anomalous dimension. At four loops the resulting anomalous 
dimension is in complete agreement with our direct computation of this quantity. Their 
proposal is compatible with integrability and with the KLOV transcendentality principle 
for the cusp anomalous dimension, and violates perturbative Berenstein-Maldacena-Nastase 
(BMN) scaling [4] starting at four loops. At five loops their proposal also matches our 
prediction of this coefficient, given in eq. (17.251) . based on using both KLV [3J] and Fade 
approximations. At higher-loop orders, their proposal corresponds to eq. (17.3 ip and appears 
to properly incorporate the numerically largest contributions. That is, it matches reason- 
ably well our approximate expression (17.301) . However, we find that successive KLV and 
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Pade approximations, based on truncations at increasingly higher orders through 13 and 11 
loops respectively, do not match the strong-coupling string results. This indicates a tension 
between the weak- and strong-coupling results, and suggests to us that further modifications 
may be necessary. The question merits further study. 

We remark that, assuming the KLOV conversion principle, our result for the leading- 
color four-loop cusp anomalous dimension also predicts a piece of the corresponding result in 
QCD. Which piece? At three loops and below, the QCD result is a polynomial in the SU{Nc) 
Casimir operators Ca and Cp, while the MSYM result is composed solely of Ca = Nc, so 
it has no subleading-color terms. The MSYM result provides one constraint on the leading- 
transcendentality parts of the coefficients of the color factors in QCD, such that, after the 
group theory Casimirs have been set to the values Cp = Ca = Nc, the leading-color terms 
are equal to the MSYM result. Starting at four loops, however, there are color factors that 
cannot be reduced to polynomials in Ca and Cp- The relevant color factors are those of 
L-loop propagator diagrams. 

In MSYM, any triangle subdiagram leads to a group-theory factor of Ca, times a lower- 
loop group-theory factor. So the question is, when do no-triangle propagator diagrams first 
appear? At three loops, there is one nonplanar no-triangle propagator diagram, but its 
color factor vanishes identically using the Jacobi identity (see e.g. ref. 52]). At four loops, 
fig. M^d^) illustrates the unique such planar graph (when the two pairs of external lines 
on the left and right sides of the graph are each replaced with single lines). There are a 
number of nonplanar graphs as well. Hence the cusp anomalous dimension in MSYM can 
now have subleading-color terms. Presumably the KLOV principle will apply with respect 
to the leading-color, planar terms, once the fermionic color factors in the QCD result are 
shifted from the fundamental to adjoint representation. But will it also apply to subleading- 
color terms? If conformal invariance is the main issue 35|, then it should. But if planarity 



is important, perhaps it will not. Of course, the question is a moot one until the cusp 
anomalous dimension in QCD, and the subleading-color part of that in MSYM, are both 
known at four loops. Similar remarks apply to the four-loop form-factor quantity Qq^^ ■ 

In order to confirm that the four-loop cusp anomalous dimension is given exactly by the 
0{a'^) term in eq. (17.251) . it would be important to evaluate it analytically. To do so, the 
four-loop integrand (14. ip would have to be evaluated analytically through 0{e~^) instead of 
only through (9(e^^) as done here. It would also be extremely interesting to evaluate the 
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four-loop integrand through 0(e°), to exphcitly check the four-loop iteration relation fl23 



for scattering amplitudes [7, 

The program MB 6^ can be used to express the coefficients at each order in the Laurent 
expansion in e as a sum of finite contour integrals. However, the number of such integrals 
increases rapidly with decreasing inverse power of e. This property, along with an increase in 
dimensionality of the integrals, makes it harder to compute the integrals purely numerically, 
requiring a substantial increase in computational resources for the same relative level of 
accuracy. For example, the order term in Af^^-* in table [TTl has a relative precision of 
2 X 10"^, 40 times larger than that of the order term, despite having had significantly 
more computing resources applied to it. Although a brute- force computation of the 1/e and 
e° terms could probably be completed with sufficient resources, it would be considerably 
lengthier than the present computation. Accordingly, we believe that additional analytic 
work would be rather desirable before proceeding further. 

The information provided in this paper offers a guide to further progress in determining 
the integrable structure of planar MSYM, as well as in studying the transition from weak 
to strong coupling. For the cusp anomalous dimension, the striking match between the first 
two coefficients of the strong-coupling expansion, as estimated by us, using our four-loop 
result as input to the KLV approximation 3J], and as obtained from string theory 33|, 
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711, provides good evidence that we have an excellent numerical understanding of 



this anomalous dimension at any coupling. The same approximation strongly suggests that 
the correct analytic forms of the four- and five-loop perturbative coefficients are the ones 
given in eq. (17.251) . These types of approximations can also be useful for checking whether a 
given ansatz for higher-order terms in the weak-coupling expansion are consistent with the 
known string theory strong-coupling results. This should help in finding the correct integral 
equation describing the MSYM cusp anomalous dimension. Another intriguing result from 
the extrapolation to strong coupling is that the next term in the expansion {0{a^^^^)) should 
be very small and may even vanish. Finally, the remarkably good numerical properties of 
the approximate formulas indicate that the transition between weak and strong coupling is 
smooth for planar A/" = 4 super- Yang-Mills theory. 
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APPENDIX A: EVALUATING FOUR-LOOP INTEGRALS BY MB REPRESEN- 
TATION 

To obtain Laurent expansions in e for our integrals we use the Mellin-Barnes (MB) 
techm^K wh,ch has bee,> successf* appl.ed numerous calculations (see. e.g.. refs. 
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591] and chapter 4 of ref. |60l]). It relies on the identity 
1 /•^+*°° r(A + z)V{-z) 



{X + YY X^+^ r(A) 27r^' ^^^'^ 

where —Re A < /3 < 0. This identity basically replaces a sum over terms raised to some 
power with a product of factors. 

It is convenient to keep arbitrary powers of propagators and numerator factors in each of 
the MB representations. In this way, one can check the MB representations by setting some 
of the indices to zero, in order to obtain a simpler integral whose value is already known. 
Also, we can obtain the values of the two non-rung-rule integrals of fig. H] from two of the 
rung-rule integrals of fig. [2] by setting some of the indices to zero. 
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We define a four-loop integral with general indices as, 

where p, r, u, v are the four independent loop, Ui the number of indices corresponding to 
the p'j propagator and numerator factors in a given graph. For each graph (x) the labels 
correspond to the propagator and numerator labels of fig. [2J In this notation a numerator 
factor is indicated by negative indices. 

Experience shows that a minimal number of MB integrations for planar diagrams is 
achieved if one introduces MB integrations loop by loop, i. e. one derives a MB representation 
for a one-loop subintegral, inserts it into a higher two-loop integral, etc. This straightforward 
strategy provides the following MB representations for the Feynman integrals corresponding 
to graphs of fig. [21 with general powers of the propagators and irreducible numerators: 



F(^)(ai,...,ai3;s,t;e) 



e 



-a—4e 



nj=2,5,7,9,ii, 12,13 r(ai)r(4 - 09,11,12,13 - 2e) 

■1- 7 rv~) 1 1 



JJ d2 ■ ( r(2 - 09,12,13 - e - 2;i,2)r(2 - 09,11,12 - e - 2:1,3) 



{27Tiy^ J_i^ "' J_i^ \sj T{aio- Z2)T{afi- z^)T{aQ- z^)T{ai- zq) 

^ r(a9 + ^i,2,3)r(a9,ii, 12,13 - 2 + e + ^i,2,3)r(2:io - z^)T{z4 - Zi) 

r(4 - 05,8,10 - 2e + 2;i,2,3)r(4 - 04,6,7 - 2e 2:4,5,6)1(4 - 01,2,3 - 2e + 2:8,9,10) 
r(2 - 05 - 08 - e + 2:1 + 2:3 - Z4 - 2:6)r(o5,9,io - 2 + e- ^1,2,3 + ^4,5,6) 

r(o3 - z9)r(oi - 2:8) 

xr(oi2 + 2;i)r(o2 + 2:7)1(2:7 - 2:io)r(2 - 05,7 - e - 2:8,10 + 2:4,5)r(2 - Oi,2 -e + Zs,io - z-j) 
xr(2 - 04,7 - e - 2:9,10 2:4,5)1(01,2,3 - 2 + e 2:7 - 2:8,9,io)r(2 - 02,3 - e + 2:9,10 - 2:7) 
xr(o7 2:8,9,io)r(2 - 05,10 - e 2:1,2 - 2:4,5) 

xr(o5 2:4 Z5 + 2:6)r(o4,6,7 - 2 e - 2:4,5,6 + 2:8,9,10) W r(-2:j) , (A3) 

i=2,3,5,6,7,8,9 



F('^)(oi,...,oi4;.,t;6) I 



X 



11^=2,4,5,6,7,9,12 r(%)r(4 - 04,5,6,7 " 2e) 
1 /•+^~ ftY'''^ r(8-0 + Oi2 -46-^9 12) 



\[dz, 



(27ri)i2 '"7-ioo \sj r(10 - o - 5e - 2:9)r(oi3 - 2:io)r(oi - 2:2) 



r(o9 + 2:9,10)1(013 - 2:10 + 2:12)1(2 - 05,6,7 - e - 2:i,2,4)r(2 - 04,5,7 - e - 2:1,3) 

r(4 - Oi,2,3 - 2e 2i,2,3)r(Oi4 - 2:4,8,ii)r(08 - 2:6)1(03 - 2:3) 
r(07 + 2:i,2,3)r(04,5,6,7 — 2 + e + 2:1,2,3,4)1^(^:9,11 — Z^)^{2 — Oi2,13,14 — C + 2:4,8,10,11 " ^12) 



X 



r(oio - 2:7)r(4 - 08,9,10 - 2e 2:5,6,7 



62 



r(a5 + Zi^4)T{2 - a2,3 - e + ^1,3 - ^5,7)r(a2 + ^5,6,7)r(ai4 - ^4,8,11 + ^12) 
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xr(a8,9,io - 2 + e - ^5,6,7 + -^9,10,11) JJ ^(-^j) , (A4) 

j=2,3,4,6,7,8,9,10,ll,12 



F^'^Hai, . . . , ai4; t; e) = xp/, ^ ^ 

1 lj=2,4,5,6,7,8,12 ^ W)'- " 04,5,6,7 " ^Cj 
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Xr(oio + ^8,10 ~ ^4,7)r(Oi3 + 2:11 — ^9)r(2 — 08,9,10 — 6 + 2:4,6,7 ~ ^8,9,lo)r(o7 + 2:1,2,3) 

11 

xr(2 - Oi2,i3,i4 - e + 2:9,10 - 2:ii)r(o9 - 2:6 + 2:8,9) Yl ^(-^j) ' (A5) 
F(^)(oi, . . . , 014; e) = r(i. \r(A n ^ 

lli=2,4,5,6,7,9,12 ^ W)^ \^ " ^4,5,6,7 " ^ej 



X 



1 r'°°-rT. r(oi2 + ^ii)r(o9 + ^i4)r(^i4-^ii) 



{2TliY^ J_,^ ■ ■ ■ J^i^ Al ^sj r(oi - 2:2)1(4 - Oi,2,3 - 2e + 2:i,2,3)r(Oi3 - ^5,10) 

r(2 — 05,6,7 — e — 2:i,2,4)r(2 — 04,5,7 — e — 2:i,3,5)r(o5 + 2:i,4,5)r(o4,5,6,7 — 2 + e + 2:1,2,3,4,5) 

r(oii — 2:8)r(o8 — 2:4,9, i2)r(8 — 01,2,3,4,5,6,7,11,12,13,14 — 4e — 2:4,6, 7,9)r(oio — 2:7,13) 
r(2 - 09,10 - e + 2:7,11,13 - 2:i4)r(2 - 02,3 - e + 2:1,3 - 2:6,8,io)r(o2 + 2:6,7,8)r(o7 + 2:1,2,3) 



X 



r(4 — 08,9,10 — 2e + 2:4,7,9,11, 12, i3)r(oi,2,3,4,5,6,7,i4 — 4 + 2e + 2:4,5,6,7,8,9,10) r(o3 — 2:3) 
xr(o8,9,io — 2 + e — 2:4,7,9,11,12,13 + 2:i4)r(2 — 01,2 — e + 2:1,2 — 2:6,7,9) 

xr(6 — 01,2,3,4,5,6,7,11,12,14 — 3e — 2:4,5,6,7,9,10,11, i3)r(6 — 01,2,3,4,5,6,7,12,13,14 — 3e — 2:4,6,7,8,9,11,12) 
xr(2 — 08,9 — e + 2:4,9,11,12 — 2:i4)r(oi,2,3,4,5,6,7,ii, 12, 13,14 — 6 + 3e + 2:4,6,7,9,11,12,13) 

Xr(2:6,9,io — 2:i)r(oi,2,3 — 2 + € — 2:1,2,3 + 2:5,7,8,9,10) 

Xr(oi,2,3,4,5,6,7,14 -4 + 26 + 2:4,5,6,7,8,9,10,11,12,13) ^{-^j), (A6) 

j=2,3,4,5,6,7,8,9,10,12, 13,14 

p4e7 (_-\ \a( _ \S-a-4:e 

F(^)(oi,...,oi5;.,i;6)= I 



nj-=2,4,5,6,7,9,ii r(%)r(4 - 04,5,6,7 - 2e)) 
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1 f^'°° r^f[dz-^^^^" r(a9 + ^i2)r(2- 05,6,7 -e-^i,2)r(a7 + zi,2, 



(27ri)i2 y_.^ "i-ioo fJl r(aio - 2ii)r(ai - 2:2)1(4 - ai,2,3 - 2e + ^1,2,3) 

r(2 - 04,5,7 - e - ^i,3,4)r(a4,5,6,7 - 2 + e + 2:i,2,3,4)r(2 - ai,2 - e + 2:1,2 - 2:5,6) 
T{as — 2:6,io)r(ai3 — 2:7)r(a3 — 2:3)r(8 — 0,1,2,3,4,5,6,7,11,12,13,14 — 4e — 2:5,6)r(ai2 — 2:4,3) 
r(a2 + 2:5,6,7)r(2 - 02,3 - e + 2:1,3 - 2:5,7,8 )r (2:5,8 - ^i)r(ai,2,3 - 2 + e - 2:1,2,3 + ^5,6,7,8) 
r(ai, 2,3,4,5, 6, 7,14 — 4 + 2e + ^4,5,6, 7,8)r(ai5 — 2:5,9)r(4 — 08,9,10,15 — 2e + 2:5,6,9,10,11) 
xr(6 — 01,2,3,4,5,6,7,12,13,14 — 3e — 2;5,6,9,io)r(a8,9,io,i5 — 2 + e — 2:5,6,9,10,11 + ^12) 

Xr(ai5 + ^12 — ^5,9)r(2 — Oil, 12,13 — £ + ^4,7,8 ~ ^9,ll)r(2 — 09,10,15 — £ + ^5,9,11 ~ ^12) 

Xr(ai,2,3,4,5,6,7,ll, 12,13,14 — 6 + 3e + 2;5,6,9,10,ll)r(2 — 08,9,15 — 6 + ^5,6,9,10 " ^12) 

12 

Xr(05 + ^i,4)r(Oi3 - Z7 + 2;9,io,ll)r(oi2 - ^4,8 + Zg) T{-Zj) , (A7) 

3=-^ 



F^^\ai, . . . ,Oi5;s,t;e) - — g 



g4e7 (_l)a(_g^8-a-4e 



X 



X 



X 



X 



n -=1 r(%)r(4 - 01,2,3,4 - 2e)r(4 - 05,6,7,8 - 2e) 

^ ^+ioo r"^Y[dz- ' -n^i^Ortas + --11,11) 



J_i^ J_i^ jj^ \sj r(oi3 - 2:10,12) r(oii - 2:14) r(oio - 2:3,13) 



14 / ■ • ■ / i i 

J —icx> J —ioc 

r(2 — 05,6,8 — e — 2:11,12,14) r(o6 + ^ii,i2,i3)r(oi3 + 2:1,2,3,4 — ^io,i2)r(oii + 2:4,5 — 2:14) 

r(oi4 - 2;7)r(o9 - 2:2,8)r(oi2 - 2:9)r(6 - 01,2,3,4,9,10,15 - 3e + 2:1,2,3,4,11,13 - 27,9,10) 

r(oi4 + 2:2,5 — 27)r(4 — 05,6,7,8,11,13,14 — 2e — 22,3,4,5,11,13 + 27,io)r(o9 — 2:2,8 + Zq) 

r(ai5 ~ 24,ii)r(6 — 05,6,7,8,11,12,13,14 — 3e + 27,9,10 — 211,13) 
r(ai5 ~ 24,11 + 26)r(2 — 02,3,4 — e — 2:7,8, io)r(2 — 09,10,15 — e + 22,3,4,8,11,13 ~ 25) 



r(ai, 2,3,4 — 2 + e — 2:1 + 27^8,9,10) 

Xr(2 — Oi,2,4 — e — 27,9)r(05,6,7,8,ll, 12, 13,14 — 4 + 2e + 2i,2,3,4,5,ll,13 — 27,9,lo) 

Xr(4 — Oi,2,3,4,9,15 — 2e + 21,2,4,11 " 26,7,9,lo)r(o4 + 27,9,io)r(o2 + 27,3) 

Xr(2 — Oil, 12, 13, 14 — e — 2i,2,4,5 + 27,9,io,12,14)r(Oi,2,3,4 — 2 + e + 27,8,9,10) 

14 

Xr(oi,2,3,4,9,10,15 - 4 + 2e - 21,2,3,4,11,13 + 26,7,9,lo) ^{-Zj) , (A8) 

where 04,5,6,7 = 04 + 05 + 06 + 07, o = X) 211,12,13 = 211 + 2:12 + 213, etc. 

Prom these general integrals, by choosing appropriate indices, we may then obtain the 
eight integrals that appear in the four-loop amplitude, 

J(a) ^ /tF(=^)(l,l,l,l,l,l,l,l,l,l,l,l,l;s,i;e), (A9) 

j(b) ^ st''F^'^\l,l,l,l,l,l,l,l,l,l,-2,l,l,l;s,t;e), (AlO) 
X^'^) = s^iF(^)(l, 1,1, 1,1, 1,1, 1,1, 1,-1, 1,1, l;s,t;e), (All) 
j(d) ^ s3tF('i)(l,l,l,l,l,l,l,l,l,l,l,l,l,-l;s,t;e), (A12) 
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j{e) ^ s2tF(^)(l,l,l,l,l,l,l,l,l,l,l,l,l,-l,-l;s,t;e), (A13) 

X(f) = lim s't 1, 1, 1,1, 1,1, 1,1, 1,1, 1,1,-1 + 77, -l;s,t;e), (A14) 

77^0 

j{d2) ^ s2tF('^)(l,l,l,l,l,l,l,l,l,0,l,l,0,0;s,t;e) , (A15) 

j(f2) ^ sVf(^)(1, 1,1,1,1,1,1,1,1, 1,1, 1,0, 0,0; s,t;e). (A16) 

For the case of integral (f), in order to obtain a valid analytic continuation to the region 
e — s> using MB 62[], it turns out that we first need to perform an analytic continuation in 
one of the indices, by an auxiliary parameter 77. After the analytic continuation we may set 



APPENDIX B: HARMONIC POLYLOGARITHMS 



We express the leading poles in the integrals and the amplitudes in terms of harmonic 



polylogarithms (HPLs) 



generalizations of ordinary polylogarithms 



92| . Here we briefly 



63( 1 for a more complete treatment. Routines for nu- 
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86| 



summarize their definitions; see ref. 
merically evaluating HPLs may be found in ref. 

The HPLs of weight n are denoted by Haj^a2...a„{x) = H{ai,a2, 
{1, 0, —1}. They are defined recursively by. 



Atfa,{t)Ha2...aAt) 



where 



a„; x), with G 



(Bl) 



= 1?^' 
H±i{x) = +ln(l+x; 



fo{x) = - 

X 



Hq{x) = Inx , 



provided that at least one of the indices Oj is non-zero. For all Oj = 0, one defines 

Hofi,...,o{x) = ^ln" x . 
nl 



(B2) 
(B3) 

(B4) 



For weight less than or equal to four, any HPL having only the parameters = and 1 
can be expressed 63|] in terms of the standard polylogarithms 92], 



U,{z) = - / -ln(l-t), 
Jo ^ 



(B5) 
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with n = 2,3, 4, and where z may take the values x, 1/(1 — x), or —x/{l — x). (For n < 4, 
polylogarithm identities imply that not all values are required.) For the four- loop results 
presented in section IVI Al through 0(e~'^), we indeed need only G {0, 1}, and weights up 
to four. In section IVIBt however, we give analytic results for lower-loop amplitudes at the 
symmetric point {s,t) = (—1,-1) through (9(e~^), in which we encounter weights up to 
six. In the Euclidean region for the planar four-point process, namely s < 0, t < 0, m > 0, 
with the identification x = —t/s, the argument x of the harmonic polylogarithms will be 
negative. 



The HPLs obey sets of identities 



631] that allow the removal of trailing zeroes from the 



string of parameters Oj. The remaining HaT^a2...a„{x) with a„ = 1 are well-behaved as x — 0; 
in fact they all vanish there. Integrals appear in the MSYM amplitudes with both arguments 
{s,t) and {t,s), so we also need a set of identities relating harmonic polylogarithms with 
argument x = —t/s to those with argument y = —s/t = 1/x. These identities are given in 
refs. [631]. A few directly relevant examples may be found in the first appendix of ref. 
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